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Chapter 1 


INTRODUCTION 

During the period June 1974 to the present , research 
relative to the understanding and alleviation of aircraft noise 
has been carried out by the investigators with grant support 
from the National Aeronautics and Space Administration. This 
report summarizes the principal results from this research. 

Among the activities during the grant period were lab- 
oratory experiments and theoretical studies on the diffraction 
of sound by surfaces with the intention of providing basic 
information relevant to the understanding of the acoustical 
implications of the engine over wing configuration. That the 
presence of the wing below the engine may partially shield 
listeners on the ground from engine noise during flyovers has 
been the topic of a number of previous reports and papers* 
and has been the subject of investigation by Hellstrom 6 , by 
von Glahn, Goodykoontz and Wagner , by Conticelli, Di Blasi 
and O'Keefe®, by Jeffery and Holbeche®, and by Sears.*® A 
principal objective is the attainment of a rational method 
for quantitatively estimating just how much noise reduction 
would be achieved by a given design. Such a method would 
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serve as a guide in the design of future EOW aircraft and 
would enable one to make quantitative comparisons of alterna- 
tive designs. 

In order to gain some quantitative insight into the nature 
of sound diffraction by wings and to provide a data base for 
the assessment of various theoretical approaches to the over- 
all problem, a series of experiments were conducted at NASA 
Langley Research Center during the summer of 1976. These were 
carried out by Allan D. Pierce and Robin Vidimos in collabora- 
tion with John S. Priesser and other NASA personnel; the 
reduction of the data was carried out under the direction of 
W. James Hadden, Jr. In Chapter 2, a summary is given of the 
nature of these experiments and of the results. 

One of the theoretical problems presented by the overall 
topic of aircraft engine noise diffraction by wings is that 
the source of the sound is not a large number of wavelengths 
away from the diffracting surface (although in cases of 
interest the listener is). Virtually all existing computa- 
tional techniques for sound diffraction by bodies are based 
on the assumption that both distances are large, so some 
analytical development was necessary to revise existing 
theories such that they would be amenable to rapid computation 
and would give quantitative insight for cases corresponding 
to the topic of wing shielding of engine noise. The details 
of this analytical study are given in Chapter 3. 
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Another topic considered during the period of the grant 
was the effect of variable ground impedance on aircraft noise 
propagation. A pertinent question is to what extent the sound 
received on the ground is characteristic of the local impedance 
near the listener and to what extent the impedance at distant 
points affects the local reception. Chapter 3, prepared by 
Dr. Hadden, gives a theory for the scattering of spherical 
waves by a rectangular area whose acoustic impedance differs 
from that of the surrounding plane. Results of experiments 
(performed during summer 1975 at NASA Langley Research Center 
by W. James Hadden, Jr., Robin A. Vidimos, and Philip Sencil) 
concerning reflection from rectangular patches are also 
described in Chapter 4. J 

A topic related to both the variable ground impedance 
problem and that of the diffraction of noise by wings is that 
of the effects of finite surface impedance on diffraction. 
Chapter 5 is comprised of a paper by the authors written 
during the grant period which summarizes the principal results 
of an analytical study concerned with this topic. 

Chapter 6 gives a theory developed during the grant 
period for the diffraction of sound from a point source by 
a thin rigid screen in the absence of ambient flow. The work 
described there is a simple extension of work reported by 
S. Candel on the plane wave diffraction problem. (See Chapter 
6 for a listing of relevant references.) Analysis given here 
shows that a simple transformation will reduce the point source 
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problem in the presence of ambient flow to one in which there 
is no flow. The solution so derived should allow some insight 
into the influence of forward motion effects on aircraft 
noise diffraction by wings. 
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Chapter 2 


LABORATORY EXPERIMENTS ON SOUND DIFFRACTION 

The experiments performed in connection with the study of 
wing-shielding of noise were divided into three parts. In the 
first experiment (Fig. 1), the obstacle used was a thin screen, 
the source was an acoustically small driver through which 
selected pure tones were projected, the source being located 
close to the barrier. Narrow-band sound pressure levels were 
measured on a circular arc far from the edge of the screen and 
also at several locations close to the screen but well inside 
its acoustic shadow. In the second experiment the previously 
described barrier and receiver configuration was used, the 
pure- tone source being replaced by a 1 inch diameter jet. The 
third experimental configuration (Fig. 2) consisted of the 
acoustic driver, a thick straight-sided barrier with a cylin- 
drical cap, and receiver and arc centered on the junction of 

the cap and the straight side of the barrier which was nearer 

to the driver. 

The source-obstacle-receiver configuration for the first 
experiment is sketched in Fig. 3. Narrow-band pressure levels 
were recorded at the microphone positions shown. Results for 
pure tone exciation of the driver at 490, 900 and 2050 Hz with 
the driver in positions 2 (level with the top of the screen) 

and 4 (9 inches below the top of the screen) are presented in 



Fig. 1 . 


Photograph of experimental apparatus for 

first experiment. ORIGINAL PACK IB 

OF POOR O’ - M.1TY 




EXPERIMENT NO. 1 


90° 



Fig. 3. Sketch of source-receiver-screen configuration 
for first experiment (see also Fig. 1) . 


Figs. 4-6. The pressure levels for microphone positions 1, 7, 
and 8, shown in Fig. 3, are presented in Table I. Although the 
pressure levels measured at a fixed distance from the edge of 
the screen show the expected trends of increased shadowing 
effect on the screen as the frequency increases and as the 
source height decreases, we strongly suspect that these data 
were affected by transmission through the plywood screen. A 
brief calculation indicates that the coincidence frequency for 
such a panel is approximately 800 Hz. Thus, the measurements 
at the lower two frequencies mentioned above may be significantly 
contaminated by sound transmission through the screen. 

The geometric arrangement for the second experiment is 
shown in Figs. 7 and 8. The one-inch diameter jet was operated 
at pressures of 2.8 and 5 psi; one-third octave band levels 
were recorded at the microphone positions indicated in Fig. 7, 
for center frequencies 500, 1000, 2000, and 4000 Hz. The 
measured 1/3-octave band levels for the reference condition 
(Fig. 7) and in the presence of the screen (Fig. 8) are com- 
pared in Figs. 9-11. It should be noted that the results for 
1000 Hz in Fig. 9 and for 4000 Hz in Fig. 10 have been shifted 
upward by 10 dB for convenience in presentation. Similarly, 
the results for 2000 Hz in Fig. 12 have been shifted downward 
by 10 dB. As in the first experiment, it is likely that 
transmission through the plywood screen is a contaminating 
artifact of the measurements in the bands centered at 500 and 
1000 Hz. The measured 1/3-octave band levels for microphone 
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Fig. 4. Measured narrow-band pressure levels 3.5 ft. from top 

of screen at 490 Hz for two souce positions (see Fig. 1) . 


P0S.2 
POS. 4 

Fig. 5. 



* WITH 10 dB DOWNWARD ADJUSTMENT 


Measured narrow-band pressure levels 3.5 ft. from top of 
screen at 900 Hz for two source positions (see Fig. 1). 


13 




’WITH 10 dB DOWNWARD ADJUSTMENT 


Fig. 6. Measured narrow-band pressure levels 3.5 ft. from top of 
screen at 2050 Hz for two source positions (see Fig. 1). 
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♦WITH 10 dB UPWARD ADJUSTMENT 


Fig. 9. Measured 1/3-octave band levels at outer 
frequencies shown for jet noise: jet 

pressure, 2.8 psi; distance from top of 
screen, 3.5 ft. (see Figs. 7 and 8). 





♦WITH 10 dB UPWARD ADJUSTMENT 


Fig. 10. Measured 1/3-octave band levels at outer 
frequencies shown for jet noise; jet 
pressure, 2.8 psi; distance from top of 
screen, 3.5 ft. (see Figs. 7 and 8). 
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Fig. 11. Measured 1/3-octave band levels at outer 
frequencies shown for jet noise: jet 

pressure, 5.0 psi; distance from top of 
screen, 3.5 ft. (see Figs. 7 and 8). 



1/3 OCTAVE BAND LEVELS (dB) 
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Fig. 12. Measured 1/3-octave band levels at outer 
frequencies shown for jet noise: jet 

pressure , 5.0 psi; distance from top of 
screen, 3.5 ft. (see Figs 7 and 8). 


Table I. Narrow-band Pressure Levels Close to the Screen In the Acoustic Shadow of a Point Source 


Frequency 490 Hz 


900 Hz 

2050 Hz 

4050 Hz 

Microphone 

Driver 


Driver 

Driver 

Driver 

Locat1on a 

Position Position 

Position 

Position 

Position 

Position 

Position 

Position 


2 

4 

2 

4 

2 

4 

2 

4 

1 

59.0 dB 

65.3 dB 

84.0 dB 

85.3 dB 

82.8 dB 

83.0 dB 

88.3 dB 

86.3 dB 

7 

69.5 

66.0 

89.8 

84.0 

82.8 

84.0 

76.5 

80.0 

8 

58.8 

56.3 

81.5 

71.0 

82.8 

81.3 

76.0 

73.0 


a Refer to Figure 1 for microphone positions. 
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positions 1, 7, 8, and 9 are presented in Table II. 

The third experiment was intended to supply information 
as to the effects of a thick barrier and a curved diffracting 
surface. The source-barrier-receiver geometry for this 
experiment is sketched in Fig. 13. The sides of the barrier 
were sheets of 1" plywood. The cap was also constructed of 
1" plywood formed so as to produce a half-cylinder with a 
radius of 12 inches. As in the first experiment, pure- tone 
excitation was applied to an acoustically small source. The 
source was located close (in terms of acoustic wavelengths) to 
one side of the obstacle. Several source heights relative to 
the highest point on the barrier were used. Narrow-band sound 
pressure levels were measured on an arc at a fixed distance from 
a point near the junction between the straight and curved por- 
tions of the barrier. Additional sound level measurements 
were made in a vertical plane in the acoustical shadow of the 
barrier at a horizontal distance of 88 inches from the source. 
The measured pressure levels for several source heights are 
presented in Tables III-V. These measurements show the expected 
increase of the shadowing effect with frequency and, in the 
main, the expected increase of the shadowing effect with 
difference between the source heights and the highest point o 
the obstacle. In some cases the variation in pressure level 
with angle is not a uniform decrease from the position almost 
directly above the source to that well inside the shadow of 
the barrier: the deviations which arise are no doubt due to 


Fig. 13. Sketch of source-barrier-receiver geometry for 
third experiment (see also Fig. 2). 





Table II. One-third Octave Band Levels Close to the Screen in the Acoustic Shadow of a 1-inch Jet 


Frequency 500 Hz 1000 Hz 2000 Hz 4000 Hz 


Hi crophone 
Location 3 

Screen 

Absent 

Screen 

Present 

Screen 

Absent 

Screen 

Present 

Screen 

Absent 

Screen 

Present 

Screen 

Absent 

Screen 

Present 

Pressure: 

2.8 psi 








1 

69.8 dB 

70.0 dB 

74.2 dB 

74.0 dB 

76.0 dB 

75.6 dB 

74.2 dB 

74.0 dB 

7 

64.5 

62.6 

68.6 

64.0 

72.0 

65.0 

73.5 

£3.6 

8 

62.8 

60.8 

65.5 

60.2 

70.0 

60.2 

72.2 

59.2 

9 

61.0 

58.2 

63.0 

56.0 

66.0 

57.8 

69.5 

56.5 

Pressure: 

5.0 psi 








1 

78.0 

77.5 

83.5 

84.0 

87.8 

88.0 

85.5 

85.0 

7 

71.0 

70.5 

77.0 

72.0 

83.0 

75.0 

85.0 

74.2 

8 

69.0 

67.5 

73.8 

68.0 

80.0 

70.5 

83.5 

70.0 

9 

57.6 

66.0 

61.2 

63.8 

66.6 

68.2 

70.8 

67.0 


3 Refer to Figs. 7 and 8 for microphone positions. 



Table III. Narrow-band Pressure Levels for Diffraction of 
Sound by a Cylindrically Capped Barrier: Source 
12" Below Highest Point on Barrier. 


Microphone 

Location* 


Frequency 


490 Hz 

900 Hz 

2050 Hz 

4050 Hz 

1 

64.5 dB 

86.8 dB 

81.5 dB 

91.5 dB 

2 

64.0 

88.0 

83.5 

90.3 

3 

59.0 

82.0 

79.8 

82.5 

4 

55.5 

80.3 

66.3 

66.8 

. 5 

52.3 

68.5 

60.5 

66.5 

7 

50.5 

74.5 

69.5 

62.0 

8 

39.5 

71.5 

58.8 

64.5 

9 

51.5 

70.0 

47.8 

70.3 


a Refer to Fig. 13 for microphone positions. 



Table IV. Narrow-band Pressure Levels for Diffraction of 
Sound by a Cylindrlcally Capped Barrier: Source 
6" Below Highest Point on Barrier. 


Microphone 

Location* 


Frequency 


490 Hz 

900 Hz 

2050 Hz 

4050 Hz 

1 

66.0 dB 

91.0 dB 

89.0 dB 

95.0 dB 

2 

63.8 

86.8 

83.3 

87.5 

3 

59.3 

84.5 

77.0 

87.5 

4 

58.0 

78.3 

69.5 

75.8 

5 

53.8 

68.3 

65.0 

70.5 

7 

49.5 

73.5 

68.0 

67.5 

8 

45.0 

75.3 

67.5 

69.8 

9 

51.3 

74.8 

65.3 

61.3 


a Refer to Fig. 13 for microphone positions. 


Table V. Narrow-band Pressure Levels for Diffraction of Sound 
by a Cylindrlcally Capped Barrier: Source at Height 
of Highest Point on Barrier. 


Microphone 

Location 9 


Frequency 


490 Hz 

900 Hz 

2050 Hz 

4050 Hz 

1 

66.3 dB 

91.3 dB 

87.3 dB 

85.8 dB 

2 

64.5 

84.5 

86.3 

89.5 

3 

61.8 

86.8 

81.5 

86.8 

4 

56.0 ’ 

75.8 

68.5 

78.8 

5 

55.0 

71.5 

68.3 

69.8 

7 

49.3 

75.3 

72.3 

74.8 

8 

47.5 

77.0 

66.3 

61.0 

9 

51.0 

71.3 

63.3 

65.0 


a Refer to Fig. 13 for microphone positions. 
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constructive interference between waves transmitted directly to 
the receiver and those reflected from the cylindrical cap. 



Chapter 3 


THEORY OF SOUND DIFFRACTION 
AROUND SCREENS AND WEDGES 
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INTRODUCTION 

Solutions corresponding to constant frequency sound dif- 
fraction by a rigid wedge or a rigid screen (a limiting case 

1 2 

of a wedge) are well known. ’ In particular, the exact so- 
lution for the case of a point source in the vicinity of such 
a wedge or screen appears in various places in the literature 
as a contour integral in the complex plane with an integrand 
of moderate complexity involving elementary transcendental 

•t A 

functions. ** This integral is not directly expressible in a 
closed form, but its value when both source and listener dis- 
tances from the edge are large compared to a wavelength can 
be expressed to a uniform asymptotic approximation in terms 

p £ n 

of Fresnel integrals ' or related functions . Expansions 
have also been derived which are appropriate to the case when 
either source or listener is close (relative to a wavelength) 

O 

to the edge. 

For those situations in which one of the distances in- 
volved is neither large nor small compared to a wavelength, 
it may be necessary to perform a numerical integration of the 
contour integral (or of other integrals which would appear in 
equivalent expressions) or to sum a large number of terms of 
the expansion appropriate to the length being small compared 

to a wavelength. Such numerical integration or summation, 
however, may be slowly convergent and may be difficult to per 
form even with the aid of a large digital computer. Although 
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direct computations of this sort have been performed by Ambaud 

g 

and Bergassoli , the method they describe, while leading to 
accurate values which agree well with their experiments, is 
intrinsically limited in application to source-listener geo- 
metries in which neither location is at an extremely large 
number of wavelengths from the edge. Further, the method is such 
that severe computational difficulties would be encountered 
were the listener arbitrarily close to the shadow zone boundary. 
While one might expect such calculations to meld with calcula- 
tions using the results of a uniform asymptotic approximation, 
the match would be evident only from a direct numerical com- 
parison. 

The present chapter is prompted by the problem of estimating 
aircraft noise shielding by wings (engine-over-wing configura- 
tion) , one of the features of which is that the sound sources 
are neither very close or very far (relative to all wavelengths 
of interest) from the wing trailing edge. Research on this 
topic should be aided by the availability of a convenient gen- 
eral purpose method for the calculation of the acoustic pressure 
(i.e., the Green’s function) at an arbitrary listener location 
caused by the presence of a unit strength point source near a 
rigid wedge or screen. Ideally, the method should be based on 
a formulation which reduces directly (without excessively intri- 
cate manipulations) to know limiting cases (i.e., source on 
edge or source and listener both far from edge). 
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Fig. 1. Geometry used to describe diffraction of sound 
waves from a point source by a wedge. 
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Such a formulation, with accompanying numerical examples, 
is presented here. Furthermore, the plots included here should 
enable one, without further need of a digital computer, to esti- 
mate the sound field and the sound reduction for the important 
limiting case when the listener is many wavelengths away from 
the edge and much further than is the source (kL>>l ,rr 0 /L 2 <<l 
in the notation explained below). Discussion is also given of 

the accuracy of approximations commonly made in acoustical 
studies. 

I. GEOMETRY AND FORMAL SOLUTION 

The geometry appropriate to the problem under consideration 
is that of a rigid wedge whose edge lies along the z-axis (Fig. 

1) in a cylindrical coordinate system (r,6,z), with the two 
faces taken as the e - 0 and e ■ 0 planes, such that the region 
exterior to the wedge extends from e * P to e * 0 (with 0 >it). 

A thin screen corresponds to 0 ■ 2w. (Here we use the same 
notation as was used in a previous paper^ by one of the authors.) 

The source of sound is a single harmonic point source (ang- 
ular frequency w, wavenumber k « w/c) located at a point (r 0 ,e o 

z 0 ) and of strength such that the acoustic pressure field p in 

i kR 

the source's immediate vicinity is given by e /R plus bounded 
terms when R, the net distance from the source, is substantially 
less than the distance of source from edge. Here a customary 
time dependent factor of e’^ ut is understood but omitted for 
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simplicity. The acoustic pressure field dependence thus cor- 
responds to a Green's function G($l£ 0 ) which satisfies the 
scalar Helmholtz equation with the customary source term 
-4ir6(j£-£ 0 ) on the right hand side. Boundary conditions cor- 
responding to the rigid wedge are that 3G/06 • 0 at 6 * 0 and 
9 ■ B, respectively. 

For present purposes, it is convenient to take the solu- 
tion to the problem just posed in the form (but in the present 

^ Q 

notation) utilized by Ambaud and Bergassoli . This, with some 
paraphrasing of notation, can be written 


G(JSL£ 0 ) - 2^ [6(c 1 )U€«-c i ) + v^)] 

i=l 


where 


t 2 - 20-|e-e o | 

<3 * e + 9 o / 

c 4 - 2 e-(e*e 0 ) 


(i) 


(2a) 

(2b) 

(2c) 

(2d) 


Here U(c) is the Heaviside unit step function. The G(c i )U(ir-c j) 
terms for i » 1,3,4 correspond to waves inferred from purely geo- 
metrical acoustical considerations, i.e., (i*l) a direct wave, 
(i»3) a wave reflected from the e » 0 face, and (i*4) a wave re- 
flected from the 0 » B face. (The term GU 2 )U(*-c 2 ) is always 



zero, since ; 2 * s always greater than it, but is included to 
preserve the symmetry of the expression.) The term G ( c ) 
represents a radially symmetric spherically spreading wave, 
generically denoted by e /R , where (arbitrary argument c) 

R ■ [r 2 + r 2 + (z-z Q ) 2 - 2rr Q cos c ] ** (3) 

This distance, for the four particular values of c listed 
above, may be interpreted as: (i»l) distance from source; 

(i=2) distance from an image-image point; (r Q , 2(e-n) + 0 Q ,z o ) 
if 0>0 Q ; (i=3) distance from image of source reflected through 
0=0 plane; and (i=4) distance from image of source reflected 
through 0 « $ plane. (While the geometrical interpretation of 
C 2 way seem irrelevant since U(ir-t 2 ) I s always zero, the inter- 
pretation is germane to the interpretation of V(c 2 ) * n the 
limiting case, termed the Fresnel number approximation, below. 

The image- image is formed either by reflecting the source 

through the 0=0 plane, then reflecting this image through 
the 0=6 plane or by carrying out the reflections in inverse 
order. The construction is indicated in Fig. 2.) In the cases 
i • 1,3,4, the presence of the Heaviside unit step functions as 
factors in the geometrical acoustics terms insures that: (i=l) 

the direct wave is zero unless the source may be "seen" by the 
listener; (i*3) there is no contribution from a wave reflected 
from the 0=0 face unless one can construct a specularly re- 
flected ray going from source to face to listener; and (i=4) 
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there should be an analogous ray reflected from the 0 = 6 face 
connecting source and listener if the corresponding geometrical 
acoustics term can contribute to the field. 

The sum of the terms V(c^) in Eq. (l) may be interpreted 
as the diffracted wave. Each may be written in a similar 
fashion as a definite integral, which, in the form taken by 
Ambaud and Bergassoli, is 


with 


V(c) 



G (ir+iw) Q (w, v, C ) dw 


Q(w,v,c) * 


(v/2) sinfv(ir-s)l 

cosh(vw) - cos[v(ir - O] 


(4) 


(5) 


the index v b*eing ir/e (v * 1/2 for the thin screen, 2/3 for a 
right angled wedge). Here G(n+iw) represents the wave function 
e ik V R being given by Eq. (3), with s replaced by n + iw, or, 
equivalently, with cos c replaced by -cosh w. The quantity ,<? 2 
is real and positive, R being understood to be the positive 
square root of R 2 , throughout the integration over w. 


II. REFORMULATION OF DIFFRACTION INTEGRAL 

Direct numerical evaluation of V(c), while possible, is 
unwieldy because of (1) the infinite limits, (2) the oscilla- 
tory nature of the integrand and the attendant slow convergence 
in many cases of interest, and (3) the fact that Q is unbounded 
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near w * 0 as e-*-ir. To avoid such difficulties we change the 
variable of integration and the path of integration. To this 
purpose, we note that Q = dij»/dw where ♦ is such that 

tan ♦ = tan[A(?)]tanh[(v/2)w] 

and where A is ( v/ 2 ) ( - 0 - t r + c) plus any multiple of tt. If we 
refine the definition of A (?) and i p such that varies from 
0 to A as w varies from 0 to », the proper choice for A is 
(given 0<t<28) 

« 

A( c ) = ( v/2) ( -6-ir+t ) + ttU(it-c) (6) 

The value of $ corresponding to its tangent as given above is 
understood to lie between and * and to have the same sign as 
A. One may note that A (c) is discontinuous at x, - n : A(c) in- 

creases from a positive value (v/2)($-ir) at c = 0 up to ir/2 at 
5 = it’, then drops abruptly to -tt/2 at c = i + and subsequently 
increases linearly, passing through 0 at ? = B+tt , up to the 
original value (v/2)(b-*) when c * 26. 

Some indication of the variation of values of the A(t^) 
[abbreviated A^ here] with the source and listener coordinates 
and e may be obtained if one considers the specific case 
■(typically of greatest interest) in which the source is on the 
far side of the wedge, 6 > & Q > it, the listener is in the 
shadow zone, 0 < e < 6 q - tt (See Figure 3). In this case all 
the A^ are negative and between -tt/ 2 and 0, the magnitudes 


FUNCTION A(Ci) 


LISTENER ANGLE 0 
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|Aj| and | | increasing with increasing e and conversely for 
|Aj| and |A 2 |. At 6 * 0, A^ » Aj and A 2 * A^; in general one 
has [AjJ > jAj| > |A,| and |AjJ > (A^J > |A 2 |. One may note 
that the line, A^ versus 0 equals -u/2 at the shadow zone boundary. 
The lines A^ and A^ cross only if e Q > (b+ti)/2 and, when they do, 
they cross at $ = 6 - e with the mutual 

value A^ = A4 = -ir/2 + (v/2)(0-ir) = -uv/2. 

If we now change the variable of integration to q = 4</A, 
then Q dw = A dq and q varies from 0 to 1. The remainder of the 
integrand can also, after some algebra, be expressed in terms of 
q rather than w. The pertinent intermediate result is 

r « [L 2 + rr 0 (Y - Y^) 2 ] 1 * (7) 

vvhere we abbreviate 

L = f(r + r Q ) 2 + (z - Zq) 2 ]* 5 ( 8 ) 

* • ) tSnUi - tan { qx l | < 9 > 

The quantity Y, and therefore the spherical wave factor, is in- 
dependent of the sign of A. Thus we may rewrite the integral in 
Eq. (4) as 


VU) » - (1/it) A(c)(e ikL /L) F v (|A|,o,e) (10) 


where 



dq 


F v ( I A | ,<*,e) 


(11a) 
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a * k rr^/L; e * rr /L 2 
0 o 

I ( q ) - (L/fl) e lk( ' ? - L) (lib) 

with L and R as given above. 

The set of arguments of F^ is readily seen from the above 
equations to be complete. The forms chosen for the parameters 
e and* a are particularly convenient in the consideration of limit- 
ing cases. From geometrical considerations, e is always less than 

1/4. The parameter a, which has the appearance of a Fresnel wave 

% 

parameter, may in principle have any value. The quantity L has 
the important geometrical interpretation of being the length of 

the shortest two segment path which goes from source to edge 
and then to listener (i.e., L is the length of a diffracted 
ray path) . 




III. THE DEFORMED CONTOUR 

The variable q is now considered as a complex variable 
and the integral over I(q) in the definition of F^ above is 
interpreted as a contour integral in the complex q plane. 

Rather than integrate directly along the real axis, we choose 

a path C which (1) goes from 0 to 1, (2) has finite length, 

(3) is such that Re(fl-L) * 0 at every point on the path, and 

(4) is such that, for nonzero a, e ' ' decreases monotonic - 

ally from 1 to 0 as q travels the path C from q * 0 to q ■ 1. 

That a path with these properties exists is supported by the 

mathematical foundations of the method of steepest descents and 
is substantiated by the construction given below. 

The evaluation of the integral along the contour C is 
facilitated by a reformulation of the function I(q), Eq. (12), 


■’WWWWW 
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The restriction Re(fl-L) * 0 along the path implies that we may 
introduce a real parameter K such that, at any point on the 
path, R is related to K by 

R = L[1 + ieK 2 ] (13) 

Here K ranges from 0 through positive values when q ranges from 
0 through successive points on the path. The relationship 
between q and K may be determined by equating the squares of 
Eqs. (7) and (13), then inserting the expression (9) for Y, 
and solving for q. In this manner one finds 

q = — — tan [tanh — tan | A ( J (14a) 

|A| 2 

with 

1/2 

sinh X * K[i/2 - eK 2 /4] (14b) 

The several ambiguities in the definitions of the square root 
and of the implied inverse trigonometric functions are 
resolved by the requirement that q vary continuously from 0 
to 1 (although not on the real axis) as K varies from 0 to «. 

To accomplish this, one defines the square root in Eq. (14b) 
to be such that its phase is between ir/4 and n/2, then defines 
X to be such that Re(X) > 0, 0 < lm(X) <_ tr/ 2 , and q to be such 
that it lies in the first quadrant. 

The computation of q R and qj for given values of K is 
generally facilitated by reducing Eqs. (14) to explicit equations 
involving only elementary functions of real variables. Such a 
reduction yields, for example 
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in which 


with 


tan(2|A|q R ) - ( 2 I A 1 ) sinh a (15a) 

cos b + cosh a cos(2|A|) 


sinh (a/v) 

. , u/ » K[(l + Q 2 ) 1 / 2 ± Q] 1 / 2 

sin (b/v) 


(15b) 


Q = \ K 2 [1 - e + \ e 2 K 4 ] 


(15c) 


The expression for tanh (2 | A | q j ) is similar to Eq. (16a): 
sinh a, cosh a and cos b should be replaced by sin b, cos b, 
and cosh a, respectively. The restrictions mentioned above 
concerning phases and btanches imply that b/v is between 0 and 
ir/2 for K < (2/e) 1 / 4 and is between tt/ 2 and it for K > (2/e) 1 / 4 . 
The restrictions further imply that 2|A|q R lies between 
0 and ir . 

Some computed plots of the deformed contour C in the com- 
plex q plane and of the corresponding variation of K along the 
contour are shown in Figs. 4 and 5. Analysis of the equations 
given above indicates that such contours always proceed from 
q * 0 obliquely upward at an angle of 45° with the real axis 
and this is confirmed by the computations. The terminal point 
q * 1, is approached from above and to the right, making an 
angle (l-v)ir with the real axis to the right of q = 1 for non- 
zero c. In the limiting case of a screen, v * h, the contour 
terminates at a right angle with the real axis. In the limit 



IMAGINARY PART 



Fig. 4. Typical deformed contours in the complex q-plane which correspond to 
paths of steepest descent for a factor in the exponential in the 
integrals described in the text. 





REAL PART 



Fig. 5. Relationship between real part of the complex variable q and the factor 
K occurring in the exponent of the integrand. (These in conjunction 
with the curves in Fig. 4 give an indication of the manner in which K 
varies along the integration contour.) 



44 


of vanishingly small e, the contour C approaches a limiting form 
which approaches q ■ 1 obliquely downward from the left, making 
an angle of v*/ 2 with the real axis. The principal modification 
of this limiting form caused by nonzero e is a small "kink” 
near q * 1 in which q R overshoots q R » 1 slightly (except for 
v * H) » the contour then bending back and approaching q « 1 
obliquely downward from the right. The quantity K always in- 
creases monotonically from 0 to • along the contour, except for 
the limiting case where |A| is identically ir/2. If |A| is 
slightly less than this upper limit, K remains virtually zero 
along the major bulk of the contour but increases rapidly to •» 
near the very end of the path. 

At this point, we may note that the reformulation of the 
diffraction integral as represented by Eqs, (10-12), with C 
taken as the integration contour, has removed all the difficul- 
ties pointed out at the beginning of this section. The limits 
of integration are now finite, the modulus of the integrand I (q) 
is bounded by 1, and the integration along C removes the problem 
of the oscillatory nature of the integrand. 

IV. LIMITING CASES 

The formulation as presented leads either directly or with 
minor mathematical manipulation to a number of important limit- 
ing expressions for the Green's function and for the various 
terms which contribute to it. 
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1. Source or listener on edge . This case is characterized 
by c * 0 and R * L for all values of q, so we have 

F v ( |A| ,0,0) - 1 (19a) 

and the total Green's function reduces to 

G(x|x 0 ) » 2vL _1 e ikL - (2 tt/ 6 )L" 1 e ikL (19b) 

where, in this instance, L is simply the distance from source to 
listener. The above pressure field, except for the limiting case 
of a thin screen (where 8 = 2ii) , is always larger than what would 
be expected were the wedge not present. The Green's function for 
source or receiver on the edge could also be derived from simple 
symmetry arguments (the field must exhibit spherical symmetry for 
source on edge, the total volume velocity of the source must be 

the same as in the absence of the wedge, but the volume velocity 
per unit solid angle increases by a factor of 4 tt / 2 6 , where 28 is the 

solid angle external to the wedge about a point on the edge) 

without the necessity of the general solution. 

2. The limi t |A| * ir/2 or { > it. In this case the ap- 
proximation R » L is valid over most of the length of the contour 
C, the contribution from portions of the contour where this ap- 
proximation does not hold becoming increasingly negligible as { A J 
becomes progressively closer to */2. Thus, we obtain 


F v (*/2,a,e) * 1 


(16a) 


so the sum of the corresponding geometrical wave G(c)U(ir-c) and 
the appropriate diffracted wave term V(c) should have the limit 

lim {G(c)U(w - c) + V(c)} - (l/2)e ikL /L (16b) 

C-*-ir 

regardless of from which side the limit is approached. Thus, the 
total field, as expected, is continuous. 

3. The uniform asymptotic limit , where krr Q /L >> 1, |A| is 
arbitrary . This corresponds to both kr and kr Q being large and 
|z - z Q | being less than or comparable to (r 2 + r 2 )* 5 . Equiva- 
lently, both source and listener are far from the ,-dge and the 
angle between the edge and the broken ray from source to edge to 
listener is not close to 0. 

In the evaluation of this asymptotic limit, it is convenient 
to regard K as the variable of integration. The derivative dq/dK 

may be evaluated by implicit differentiation of Eqs. (15b) such 
that dq/dK is a function of a and b times the derivative 
d(a+ib)/dK. Since krro/L is large we may expect the dominant 
contribution to the integral to come from small values of K. 

However in the limit K -*■ 0, dq/dK is inversely proportional to 
cos(2|A|) and is singular when |A| ♦ ir/2. To cover this contingency 
one expands the denominator in the function just mentioned to the 
next order nonvanishing term (which turns out to be second order) 
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in K. The remainder of the factors (except for the exponential) 
are approximated by their limits as K + 0, In particular, one may 
note from Eqs. (15) that d(a+ib)/dK is just v(l+i) in this limit. 

The variable of integration is next changed to u * a 1 ^ 2 K, then 

the resulting integral is recognized as a constant times the integral 



* f ( X) - i g(X) (X > 0) 


(17) 


where 

X * [4ot/ir J* 4 ( 1/v ) cos ( | A | ) (18) 

Here F(X) and g(X) are the auxiliary Fresnel functions discussed 
in 9 previous paper^ by one of the authors and which are tabu- 
lated on pages 323-324 of the NBS Handbook of Mathematical Func - 
tions The mathematical manipulations as outlined above then 
lead to the expression 


F v * (*//?) e W4 [(sin|A|)/|A|][f(X) - i g(X)] 


(19) 
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for krr 0 /L >> 1. One may note that, although the coefficient of 
cos(|A|) in Eq. (18) is presumed large, it cannot necessarily be 
assumed that X is large since cos(|A|) would be very small were 
| A | close to ir/2. 

In the limit of large X, the quantity f - ig approaches 
1/ ( wX) and thus F y decreases asymptotically as the inverse square 
root of a for nonzero value of cos(|A|). When |Aj approaches ir/2, 
both f(X) and g ( X ) approach the value 1/2, the limiting values for 
X -*■ 0. In this limit goes to 1, just as indicated by Eq. (16a). 
It should also be noted that in this approximation F^ is indepen- 
dent of the parameter e for a fixed value of a, 

4. The Fresnel number approximation A 2 if, in addi t :-n to 
krr Q /L >> 1, it is true that cos(|A|) is substantially less than 
v,the parameter X in Eq. (18) may be interpreted as X 5 (2N)* 5 
where N is a Fresnel number given by 

N - (L - R a )/(X/2) C20) 

which represents the excess of the diffracted path length L be- 
yond some direct path length R A in units of half wavelengths. 

The appropriate identification of R A is 

R a - [r 2 ♦ r 2 ♦ (z~ Z Q ) 2 - 2rr 0 cos(B v )]’ s (21) 

with B v ( |A| ) taken as 
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B y (|A|) =■ + {ir - (2/v)(ir/2 - |Aj)} + 2mr (22) 

with n being an integer (0, positive, or negative) and with any 
choice of the two signs. With the purpose of giving a meaning- 
ful geometrical interpretation of B v , one may show with some 
effort that it is possible to choose the sign and the integer 
n such that 


0 

0 

u * 

l« ' 9j) 

(23a) 

0 + 2(6-n) 

0 

u - 

2 e - 1 e - e |, e > e ) 
1 0 1 9 0 

(2 3b) 

0 0 - 2 ( 8 -it ) 

u = 

28- ! 0 - 0 1 , 0 >0) 
' 0 1 f 0 

(2 3c) 

-0 

(c = 

0 + 0 ) 

(23d) 

0 

0 


26 - e o 

U = 

28 - 0 - 0 ) 
0 

(2 3e) 


Thus, with reference to the discussion following Eq. (3), is 
the direct distance of listener from (i * 1) the source; (i = 2) 
the image of the image; (i - 3) the image formed by reflection 
through the 0 « 0 plane; or (i * 4) the image formed oy reflec- 
tion through 0 * 0 plane. 

That X is approximately ( 2N) 15 where N is as defined above 
in the limits cos(|A|) << v, a, follows from the general expres- 
ion (22), from the (consistent) approximation sin[ ( 1/v ) (ir/2 - |A|)] 
* ( 1/v )cos ( ( A| ) , from the fact that e is always less or equal to 
1/4, from the definition (8) of L, and from an appropriate binomial 
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expansion of 

When the Fresnel number approximation is valid, | A 1 should 
be close to w/2, so it is consistent to approximate the 
sin( | A| )/ | A | factor in Eq. (19) by 2/w and the resulting expres- 
sion for F v becomes 

F y * (/I) e ilr/4 {f ( [2N] 5 *) - i g( f 2N1 ^ ) } (24) 

This represents a considerable simplification in that the right 
side depends on one and only one parameter N of relatively simple 
geometrical interpretation. There is no explicit v, | A | , a, or 
e dependence, other than the :uanner in which these enter into the 
determination of N. The expression above also has the virtue of 
never giving a magnitude of F v greater than 1. 

The corresponding expression for V(c) in the Fresnel number 
approximation may be obtained from Eq. (10) with A(c) replaced 
by ( tt/ 2 ) sin ( w - ?). This is in accordance with Eq. (6) and the 
fact that | A ( should be close to w/2. Consequently, Eq. (24) 

i IcT 

should be multiplied by sin(e - w)(2L) -1 e to obtain V(c). 

5 . The case when kL is^ large but krr Q /L i£ finite or t « 1 , a 
finite . The two statements are equivalent since kL -*■ • with krr Q /L 
fixed implies rr Q /L 2 ■* 0. This limiting case is of interest in 
those problems where the source is at finite or small distance re- 
lative to a wavelength from the edge but the listener is at a large 
number of wavelengths from the edge, much further than is the 
source. Conversely, because the solution conforms to reciprocity 
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(interchange of source and listener), the corresponding limiting 
solution corresponds to the pressure field in the vicinity of the 
edge when the source is a large distance away. In this reciprocal 
problem the incident wave near the edge is very nearly planar, so 
the limit can be obtained from the solution of the related problem 
of plane waves incident on a rigid wedge. The limiting case, 
source near edge, listener far from edge, is of principle interest 
in aircraft noise problems where the source is in the vicinity of 
r. wing but the listener is on the ground at a large distance away. 

The limiting value of the diffraction integral as 
rr Q /L 2 -+■ 0 may be simply denoted as F (|A|,o,0). The limit exists 
and may be readily obtained from the formulation given in the 
previous section by (1) replacing the factor L /H in the integrand 
by 1 and (2) setting e = 0 inEqs. (14) and (15). This 
yields sin(b/v) * tanh(a/v) and Eq. (15b) gives K = sinh 2 (a/v)/ 

— F 2 

cosh(a/v). The integrand I(q) reduces toe a along the contour C. 

The value of the integral F (|A|,a,0) for | A | * n/2, or for 
a =0, or for a >> 1 may be inferred from the cases 1-5 discussed 
above. Thus F^ is 1 for |A| = n/2 or for o ■ 0 and is given by 
Eq.(19) for a >> 1. Also, the Fresnel number approximation, Eq. 

(24), should be applicable in the double limit o >> 1 and cos A << v, 
the appropriate identification for the Fresnel number N in the 
limit e -► 0 being 


N - 4[rr o /(U )]cos 2 (B v /2) 


(25) 
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As regards the behavior of F v (|A|,a,0) for a << 1, one can 
derive an expansion of the contour integral in honinteger 
powers of a, the starting point being 

F v (|A|,o,0) = 1 -jf (1 - e- ctK )(dq/da)da (26) 

In view of the restriction krr Q /L << 1, the first factor in 
integrand above is small unless a is relatively large. Thus, 
if we seek just the leading term and anticipate that this, for 
sufficiently small values of the expansion parameter, is larger 
than any given constant times this parameter, it is sufficient 
to adopt the approximations K 2 = (l/2)e a ^ v , dq/da * | A | _1 sin 
(2 1 A| )e iv11 / 2 e' a , i.e. asymptotic limits for e = 0, a large. 

Then the variable of integration may be changed to u * (1/2) a e a ^ v 
such that (dq/da )da is a product of u- independent factors and 
u" v_1 du, one of these factors being [a/2] v . The lower limit on 
the u integration becomes a/2, but, providing v is not very 
close to 1 (i.e., we here exclude the case of highly obtuse 
wedges) , this can be approximated by 0 insofar as we are only 
interested here in the lowest order (which is lower than first 
order) term in a. In this manner, one obtains 

F v (|A|,a,0) * 1 - |A|- 1 sin(2|A|)e _ivn/2 [a/2] v r(l-v) (27) 

Here we recognize (after integration by parts) that the inte- 
gral over u of v(l-e~ u )u’ v * 1 is the gamma function with ar- 
gument 1-v. 
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The fact that v is less than 1 implies that the magnitude 
of F v decreases sharply from 1 (the derivative of its magnitude 
with respect to the expansion parameter is negative and becomes 
singular when the parameter approaches zero) when a increases 
from zero. As discussed subsequently below, this implies that 
a modest amount of sound reduction in the shadow zone is 
achieved even when the source is only a slight distance from 
the edge. 

In this same limit of rr Q /L 2 -*■ 0, krr Q /L << 1, the total 
Green's function (found by inserting the above into Eq. 1) 
becomes 


G(x|x Q ) « (2w/e)L’ 1 e ikL {1 + 2e' iv7r/2 [1/r ( 1 + v ) ][krr o / (2L) ] v 

(28) 

cos(ve)cos(v0 ) ) 
o 

where we make use of the identity 

sin(vu)r (l-v) * vn/r(l+v) 

The above approximate Green's function is consistent with a 

g 

more general expansion given by Tuzhilin. One may note that, 
if the listener is in the shadow zone, cos(ve) and cos(ve o ) 
have opposite signs, so the second term in Eq. (28) would de- 
crease the magnitude of the Green's function in such cases 
(as should be expected) from that represented by just the first 
term. The phase of the Green's function is predicted to be 
greater than kL. (The formulation in general requires the 
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phase in the shadow zone to lie between kL and kL + ir/4.) 

6. The case of a thin screen (v = 1/2) for e -*• 0 with 
a finite . For the most part, it is conceptually simpler to 
consider each V(c^) in Eq. (1) as being calculated individu- 
ally, the sum being found subsequently. Although these occur 
in pairs, V(?) and V(2g-e), there appears in general to be no 
major analytical simplification obtained by considering such 
a pair as a unit. An important exception is the case of the 
thin screen (v = 1/2). The fact that some simplification 
should be possible in this limit should be evident from the 
fact that the geometry of source, images, and image- image in 
this limit is degenerate: the source and image-image coincide 

and the locations of the two images coincide. The analytical 
simplification is of minor computational advantage except in 
the limit e -*• 0. The simplification which results in this limit 
(which, as pointed out above, is equivalent to the problem of 
diffraction of plane waves by a thin screen) is that the 
Green’s function and each of its two constituent pairs, 

V(?|) + V(2g-Cj) and Vfej) + V(2g-?j), can be expressed rather 
simply in terms of Fresnel integrals. (Given the incident 
plane wave interpretation of this limit, this is a well known 
result .) 

The manner in which the result may be obtained from the 
formulation presented here is first to change the integration 
over q to one over a. Then the sum V(c) + V(2g-c), with V(c) as 
given by Eqs. (10-12), with the q integration along the contour 


SB 
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C, may be grouped as a single integral over a from 0 to » 
which involves a factor 


A(c ) dq( | A c ) | »a)/da + A(2&-c) dq({A 2B-c)|,a)/da 


One should note that q, considered as a function of A and a, 
will in general have different values if } A f is taken as 
|A(c)| or |A(28'c)|* Evaluating this expression for v * 1/2, 

$ = 2ir, e = 0, such that sin[2| A(c) | ] * |cos(c/2) | , cos [2 |A(c) | ] * 
-sin(c/2), tan b ■ tanh a, K 2 » sinh(2a)tanh(2a) , etc., it 
eventuates, after some lengthy algebra and application 

9 

of various trigonometric identities, that this can be expressed 
rather simply as a function of K and cos(c/2) times the der- 
ivative dK /da with no explicit dependence on a. Consequently, 
the variable of integration can readily be changed to u * a 1 / 2 K. 
Once this is done, the integral appears in the form of a constant 
times the diffraction integral A D (X) of Eq. (17) with the ap- 
propriate identification for X being 

X * (,4a/ ir]**! cos( t/2) | <- 2g 

In this manner, we obtain 


A(c)F lj {|A(c)|,o,0) + A(28 -c)F, s (|A(28-c)|,o, 0) 
-sign{cos(c/2)}(n/2 ,s )e iir/4 [f(X)-i g(X)] 


( 30 ) 


with X as given above. The corresponding expression for 
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V(c) + V(2B-c) is just -(l/w)L _1 e lkL times Eq. (30). The 
total Green's function may then easily be written down from 
Eq. (1). In the case where the listener is in the shadow zone 
(diffracted field only), cos(c/2) is negative both for 
C - I e ~ e 0 1 a nd for c ■ 9 +e 0 » so the field is 


G(x|x 0 ) - 2" J *L~ 1 e i ^ L e ilT ^ 4 ([f(X) - i g(X)) t . | | 

♦ IfOO - i *(x)] ? . 6 + 0o > 

in which the indicated values of c are to be used in Eq. (29) 

J 

to compute the variable X. 


(31) 


V. NUMERICAL INTEGRATION SCHEME 

We return now to the general problem of determining the 
integral F y . The integral over l(q) along the curve C can be 
symbolically written 

F v * f l(K,e ,a)dq (32) 

•T 


where 

I (K,e ,o) - (1 ♦ ieK 2 ) " 1 e' aK2 (33) 

The quantity K is that given implicitly by Eqs. (15) and may be 
considered a monotonically increasing real function of distance 
along the contour C. 
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The prototype integration scheme suggested is one in which: 
(1) the variable of integration is first changed to K; (2) the 
domain of K integration is broken into N + 1 intervals (0, K^) , 
(Kj, K 2 ) » (Kjj, «) where N >. 1; (generally one takes 

N = 1) and (3) the integration over the first N intervals is 
transformed through an "integration by parts". Thus one has 

N 

n»l 


r 


J(K, e, o , | A | ) dK 


(34) 


K 


n-1 


+ ^ ^ * 6 * 0 ) Q ( e ) + 



(I) (dq/ dK) dK 


where 


J (K) - -2 I (K) q(K)K 

We also use the fact that q(K) * 0 if K * 0. 

One may note that the real and imaginary parts of the 
function J (K) are bounded and continuously differentiable and 
that these component parts are certainly not oscillatory. 

Thus, one may expect that the first N integrals of the above will 
be amenable to any numerical integration scheme which, while 
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utilizing values of the integrand at only a relatively limited 

number of points (less than, say, 10), achieves a high accuracy 

because of the "smoothness" of the integrand. Possible integration 

formulas (Chebyshev's equal weight, Gauss’s, or Lobatto’s, for 

example) are summarized in particular in Sec. 254 of the 

Handbook of Mathematical Functions 1 4 . (Our experience has been, 

in the present context, that 10 point Lobatto integration 

invariably gives at least eight digit accuracy.) 

As regards the integral from K N to », the qualitity 

1 1 (Kjsj) I 1 1 ~q (Kjsj) I may for most practical purposes be considered 

as an upper bound to its magnitude. It may be presumed that one 

has chosen K N sufficiently large, either that the magnitude of the 

integral is definitely negligible within the desired computational 

accuracy or else that the e-“ KZ factor in the inte * rand do " iMtes 
its decay. In the former case the last term is discarded while 

in the latter case it is evaluated by (1) integrating by parts 

and (2) performing the integration over the resulting expression, 

which has the form (representing the sum of the last two terms in 

Eq. (34). 

0,1(2 L(K) dK 

(with an obvious identification for L(K)) by Hermite integration. 14 
(Our experience is that an 8 point scheme is more than adequate) . 

The choice for the K^, ...,K^ as well as the parameter N 
should not be too critical. One could compare answers obtained 
with different choices of these parameters in order to assess 
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whether or not some desired accuracy has been obtained. One 
could, for example, simply take N ■ 1 and * 1/a, unless 
a were extremely small compared to unity. (We have at present 

a somewhat elaborate scheme for chosing these parameters, but 
the details seem too arbitrary and unimportant to warrant their 
inclusion here.) 

Computation time for a single value of may be considered 
as roughly directly proportional to the number of times which the 
function q(K) must be computed from Eq. (15) (which is a straight- 
forward evaluation requiring trigonometric and inverse trigonometric 
functions). This number is typically just 18 with the scheme 
as outlined above so the computation time should be of minor 
consequence, given the availability of a modern high speed 
digital computer. 

Some sample calculations are presented in Figs. 6 and 7. 
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RELATIVE DISTANCE PARAMETER krr./L 

Fig. 6. Magnitude of diffraction function F v vs. relative distance parameter 

krr Q /L for the case v s 2/3 and for the limit rr 0 /L^ -*•(). The curves 

are shown for several values of the angular factor A. The results 
indicate that the uniform asymptotic approximation described in the 
text is valid for krr /L > 1. 


O' 
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PHASE OF Fv (DEGREES) 


ik 



Fig. 7. Phase of the diffraction function vs. relative distance parameter 

krr 0 /L for the case v * 2/3 and for the limit rr 0 /L^ -*•(), for various 

values of A. The figure indicates that the A=0 curve is a good 
approximation except for A*ir/2. 
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Chapter 4 


SCATTERING OF SPHERICAL WAVES 
BY RECTANGULAR PATCHES 


The body of this chapter consists of a copy of a 
paper prepared for submission to the Journal of Sound 
and Vibration by W. James Hadden, Jr., Robin A. Vidimos 
and Philip M. Sencil. [The experiments described in 
the paper were performed in an anechoic chamber at 
NASA Langley Research Center (Fig. i).] 
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Abstract 


A theory is presented for the scattering of spherical waves 
by a rectangular area whose acoustic impedance differs from that 
of the surrounding plane. This theory extends previous analyses 
to include diffraction effects explicitly. Results of experiments 
concerning reflection from rectangular patches are also reported. 
Agreement between these results and predicted values is not 
uniformly good, although improvements could be achieved through 
alterations in the measurement procedure. 
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INTRODUCTION 

The present paper is motivated by an interest in the effects 

of acoustical characteristics of the ground on sound originating 

in low-flying aircraft. As part of this study, analytical and 

laboratory investigations have been performed on the reflection of 

sound by plane surfaces of known acoustic impedance [1*| . In 

* 

analyses of the reflection of spherical waves by i plane surface 
on which a local -reaction impedance boundary condition is imposed, 
it is customary to employ the method of steepest descents in order 
to obtain an approximation for the reflected pressure [2,3]. The 
use of this approximation can be interpreted in terms of geometrical 
acoustics as neglecting the effect of waves scattered from regions 
of the surface outside a neighborhood of the shortest reflected 
ray path from the source to the receiver. The investigation with 
which the present paper is concerned sought to determine the size 
of the effective area near the vertex of the reflected ray. This 
information could be used in developing a simplified technique for 
predicting the received sound for moving sources near the surface. 

In the interest of simplicity, experimental measurements were 
made in an anechoic chamber of sound pressure levels above rectan- 
gular patches of various areas. Pure tones were used to excite a 
small source. Sound pressure level measurements were made in the 
direction of the presumed reflected ray path. These experiments 
are described more fully in Section IV. In conjunction with the 
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experimental work, a theoretical investigation of scattering by 
rectangular areas was undertaken in which diffraction effects due 
to the finite size of patches were included. This analysis is 
discussed in Sections I-III. 


I. THEORETICAL EXPOSITION 


The analytical development is roughly parallel to that of 
Morse and Ingard for plane wave incidence [4]. The surface z - 0 
contains a rectangular patch with point impedance pcn^; outside 
the patch the normalized impedance is taken as n. The geometry 
is illustrated in Figure 1: A point source is located at ( r s »° s »4 s )> 
the receiver coordinates are (r,9,<j>). 

The received pressure may be expressed, employing Green’s 
theorem, as 



in which the Green’s function G(r|r 0 ) is approximated by terms 
representing a source point r Q = (r 0 *® 0 »^ 0 ) and a single image 
point r^ = ( r 0 » 1T ' e 0 »4> 0 ) with the image source strength (a modified 
plane-wave reflection coefficient) chosen such that the condition 
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in 3G(f|r ) 

G(r ’V-T» ' 


( 2 ) 


is satisfied to a better degree of approximation than could be 
obtained by using the plane-wave reflection coefficient. The 
approximation to the Green's function is 


G(r|r) 


ik|r - r 0 | 
e * 

4ir lr ■ r 0 l 


ik lr - r o! 

e 

4ir|r - r^| 


(3a) 


n B' cose* - 1 

R' = 

n B' cose^ + 1 


(3b) 


B' 


i 

k lr - r 0 l 


(3c) 


where is the azimuthal angle between the source- to- receiver 
point line and a line parallel to the z axis, and the inclusion 
of the factor B' represents an attempt to account for the curvature 
of the wavefront. 

The pressure terms in the integrand of equation (1) are 
approximated in a similar fashion as a combination of waves 
incident from a point source at r and an image source at r' 
below a plane characterized by the normalized impedance n^. The 
appropriate form for this approximation for the pressure may be 


inferred readily from equations' (3) with suitable modifications 
of parameters. Thus the "direct" pressure term in equation (1) 
is taken as 
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ik|r - r | 


W2 


Air, 


ik| r - r sl\ 
lr • r sl / 


(4) 


in which R g has the form of equation (3b) with 

B'+B ■ 1 + i/kr_ and e' •+• e. = cos’^z./r,.) . The scattered 

5 5 OS S S 

pressure term may be written as 




ik(|r - 
e 






B' B^ cose* cos9^ (n - n A ) 

(1 + nB' cose')(l + n A B^ cose^) 


(5) 


In order to obtain a closed- form expression for the pressure at 
some distance from the scattering area it is expedient to expand 
the factors in equation (5) which involve the distances |r - r Q | 
and |r Q - r $ | as power series in x Q and y Q . The expansions of 
such factors multiplying the exponential in equation (5) may be 
truncated so as to yield a desired accuracy which depends on 
ratios such as L/r and W/r. However, in the exponent the 



71 


criteria! governing truncation of the expansion involves the 

2 

Fresnel wave parameters, which have the form r/kL . 

Retaining second-order terms in x Q , y Q in the exponent in 
equation (5) yields approximations for the scattered pressure in 
which diffraction effects are readily discernible. In addition, 
tjiis treatment allows one to investigate the transition from the 
Fraunhofer diffraction regime (large Fresnel parameter - equivalent 
to the Morse- Ingard treatment [4]) to the Fresnel diffraction 
(small Fresnel parameter) range and beyond to the ray theory 
limit. An outline of the present expansion of equation (5) is 
given in Appendix A. The scattered pressure is approximated by 


Psc^ s 


kLW e 


ik(r + rj 


4irrr_ 


P S c I ( a i» a 2* e l»®2 ,Y ^ 


( 6 ) 


with the abbreviations 


P 

sc 



B B cose cose 

3 3 

(1 + n B cose)(l + n A B s cose g ) 


(7) 


and 


I 



1 -i (o 1 X-0 1 X 2 )r 1 -i[(a 2 +YX)Y - 0 2 Y 2 ] 

dX e Id Ye 

•'-1 


X(1 + MX + NY + QX 2 + KY 2 + SXY) (n - n A ) (8) 
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and, finally 

“l ■ < s “ 9 IS ♦ * sin9 s IS |(w) w 

2 x 7 

2 

6, - [r (1 - sin 2 e cos 2 ♦) 

1 

* r(1 ' sin \ ISz (9b) 

2 2 kLW 

Y = (r_ sin 0 sin2$ + r sin e_ sin2<j>_) (9c) 

5 s s 8 it s 

The parameters and «2 involve projections of the scattering 

area's dimensions (normalized by wavelength) on the lines from 
source and receiver for the center of the area. The parameters 
pp &2 y are similarly projected inverses of Fresnel wave 
parameters. These parameters characterize the diffraction effects 
in the approximation for the scattered pressure. The coefficients 
M, N, Q, R and S in equation (8) , in addition to providing 
correction terms depending on the size of the scattering area 
relative to source and receiver distances from the patch, are 
functions of the other geometrical and impedance parameters. The 
coefficients M and N are linearly dependent on quantities such as 
sine, sin<(> and L/r or W/r. Q, R and S are quadratic in these 
quantities. Explicit expression for these coefficients are given 
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in Appendix A. 

II. PATCH WITH CONSTANT IMPEDANCE 

For cases in which the impedance of the scattering area is 
constant, the integrals in equation (8) could be evaluated by 
completion of squares in the exponents followed by application 
of standard integration formula, but for one complication - the 
inner integral (e.g., the integration with respect to Y in 
equation (8) results in several terms involving Fresnel integrals 
[5] whose arguments have the form, in this case. 

The presence of the second integration variable precludes exact 
analytical evaluation of the remaining integration. However, 
reference to equations (9) indicates that the X-dependent and 
unity terms in the arguments are of order (L/r) compared to the 
«2 terms. In addition it can be seen that both 02 and y vanish 
in the important case of specular reflection (# * e , $ ■ 0, 

$ s * m). For these reasons , and in view of the behavior of the 
Fresnel integrals in the small- and large -argument limits [5], 
it seems a reasonable approximation to neglect the X-dependent 
terms but to retain the unity terms in the arguments exemplified 
by equation (10). 
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If this approximation is accepted and the resulting expression 
simplified by neglecting terms which are of order (kr)" 1 , (L 2 /r 2 ) 
or smaller, the integral in equation (8) may be approximated as 


I 


7T 

2(6 1 6 2 ) 1 '' 2 



AF(02 ^2/282) AjAF(u 82 ,voij/2 


* - 1Bl (j** . e i4 2- 

(2wp0^)l/2 \ 



with 


u « 1 -y 2 /48 1 0 2 , v = 1 + a 2 Y/2a 1 0 2 

(vo^) 2 al 

- + — 1 

4u8l 48 2 


(12a) 


(12b) 


•» • ->( 3 ; •■)’ 


(12c) 
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(12d) 


and 


AFfa,b) = F f [a 1/2 (b + 1)] - F £ [a 1/2 (b - 1)] (13) 


in which we have employed the abbreviation F £ (x) = C(x) + iS(x) , 

C and S being the well-known Fresnel integrals [5] . The coefficients 
of the several terms in equation (11) are 


B 


2 



(14a) 


B 


1 


M + 


Ny + S«2 

2e 2 


+ 



(14b) 


2 

Not 7 Rou va, 

A = 1 + — - + —j* —L. B 

2 ($ 2 4e£ 2y6 1 


+ 



(14c) 


Aj^ contains terms of order unity. The terms in equation (11) 

- 1/2 

which involve Bj and B 2 are of order (kr) . General 
expressions for M, N, Q, R and S are given in equations (A8). 
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One may check that equation (11) reduces to an extension of the 
result reported by Morse and Ingard [4] by noting that in the limit 
as 8^, P 2 and y become very small the function AF(a,b) [equation 
(13) ] , with argunents such as those in equation (11) , may be 
approximated [6] as 


AF(a,b) 




sin(2ab) 


(15) 


Upon substituting this egression in equation (11) , the first 

term reduces to a form similar to equation (8.3.5) of reference 

4. The second term in equation (11) vanishes in this limit, while 

1/2 

the third term is of order (e ) and hence negligible. 


III. TWO LIMITING CASES 

Although considerable simplification in the above expressions 
for the scattered pressure may be achieved in several interesting 
special geometrical configurations - forward scattering (<p = <f> s - ir) 
and specular reflection (<|> .« <p s - n and e = e s ) - only two special 
or limiting cases will be considered in detail here for brevity. 
The first, which is relevant to the experiments reported in 
Section IV, concerns reflection in the special case in which the 
source and receiver are in the plane bisecting perpendicularly 
the scattering area, i.e., 4 S * t,. [The case in which the 
source and receiver are in a plane parallel to the x - z plane 
of Figure 1 can be treated by an obvious modification of the 


limits of integration in equation (8) . ] The second case for which 
a compact expression for the scattered pressure can be obtained 
concerns scattering by a strip (taken here as lying along the y-axis 
of Figure 1). 

An explicit expression for the scattered pressure can be 
readily obtained in the special case of specular reflection with 
4> s = n: 

^ik(r+r g ) g ^ co - sQ ^ (n-n A ) 
p spec p o (r+r s ) (1+nB cose s ) (l + o A B g cose s ) 

^^4F f (6} /2 )F f (6^ /2 ) (16) 

where the factors B and B g are defined after equation (4) . In 
this case the parameters Bj and $2 become 


8 


1 

2 


k (r + V L 2 cos 2 e | 

8rr g W 2 ) 


( 17 ) 


The reduction of equation (16) to the form obtained by Leizer [7] 
for a rigid rectangle is readily apparent if one takes the limit of 
equation (16) as the normalized impedance n A becomes very large. 

An expression for the scattering by a strip of width L may be 
obtained from equations (11) -(14) by considering the limit as 
,a 2 » & 2 » and Y become very large. It is also convenient to t3ke 
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advantage of the /-translational invariance of the geometry by 
setting $ s * w. In the case of scattering in the specular plane, 
the scattered pressure term reduces to 

i [k(r+r )+1/4tt J 
V)C .strip '' 2rr 5 (8^5 X?r 




in which the parameter (<-, of equation (9b) has been modified to 

*i 


Jk 

2rr 


s 


(r ♦ r g ) 


(19) 


to produce a form consistent with the direct computation from 
equation (A4) seq . with the y -limits set to infinity. The 
coefficients M and Q in tins case are: 




(2 ♦ nB cose) 
(J * r(B COSO) 


sine! 1 

♦ I same (20a) 

J Js 
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L 2 sin 2 e (2+n B cose) (1-3 sin 2 e) 

— j - x/2 

4r 4 L(l + n B cose) z (1+n B cose) 


♦ 


l} l sine sine g 

— y [same] s + 

4r g rr s (l+n B cose)(l+n A B g cose s ) 


[2+ nB cose+n A B s cose s +l/2nn A BB s cose cose s J 

For e» e $ (specular reflection), equation' (18) reduces to 

i[k(r+r s ) + l/4it] 

fWl,strip * 5 2 p sc F f(»! /2 ) < 21 > 

(r+r g ) cose s 


with gj given by equation (17). 

IV. EXPERIMENTS ON REFLECTION 

In the experimental phase of this investigation, measurements 
of sound pressure level were made in the specular reflection 
direction above rectangular scattering areas composed, in one 
instance, from 4' x 8* (1.22m x 2.44m) sheets of 3/4 -inch (0.019 m) 
plywood laid on the floor of the Anechoic Noise Facility at the 
NASA Langley Research Center; in a second set of measurements the 



plywood was overlaid with one-inch (0.025 m) glass-fiber panels. 

In each case, pure tones were projected from a source small 
compared with the acoustic wavelength. The source and receiver 
(a 1/2- inch microphone) were arranged so that the specular plane 
bisected the scattering area. Incidence angles of 70° and 80° 
were used. Normal impedances of samples of the plywood and 
glass- fiber plus plywood were obtained from impedance tube measure- 
ments. 

The measured iirpedances were employed in computations based 
on equation (16) ; the "background" specific impedance was assumed 
as unity. The measured impedances for two selected frequencies are 
presented in Table I, Comparisons of the experimentally obtained 
sound pressure levels with those computed from equations (1), (3) 
and (16) are presented in Tables II- IX. Because the primary 
interest in this study was the variation of the reflected sound 
with size of the scattering area, all measurements have been 
normalized to the experimental result for the largest rectangle. 

As may be seen from Tables I I- IX, the agreement between 
experimental and theoretical results is by no means uniformly 
good. Two possible causes of the discrepancies are suggested: 
First, the assumption that the impedance of the grill -work floor 
which surrounded the scattering areas can be taken as that of air 
is suspect. Second, there is the possibility of a distributed- 
reaction effect in the measurements. The former question could be 
resolved by further measurements of sound pressures above the 
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bare floor of the anechoic chamber. The second possible problem 
could be rectified by the inclusion of a distributed impedance in 
the development following equations (3) . 


V. CONCLUSION 

A theory has been presented for the scattering of sound by 
rectangular patches characterized by uniform (local-reaction) 
acoustic impedances . The theory explicitly includes diffraction 
effects absent from previous analyses. Comparison between this 
theory and a set of laboratory experiments reveals discrepancies 
which may be reduced by changes in the measurement procedure or in 
the analytical model. 
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RECEIVER 


Sketch of the source-patch-receiver configuration used 
in the analysis of scattering of sound by a rectangular 
patch in which the acoustic impedance differs from that 
in the rest of the plane including the patch. 




Table I. Measured Specific Acoustic Impedances of Scattering Areas 


Frequency 

Plywood 

Glass-Fiber over Plywood 

1600 Hz 

41.89 + i56.69 

0.08 - il.03 

3200 Hz 

2.33 + i22.62 

0.10 + i0.17 

Table II. 

Relative SPL above Rectangular Areas of Plywood with 
Receiver Distance 2.7 m and Incidence Angle 80°. 


Scatterer 

Dimensions 

(m) 

f - 

1600 Hz 

f = 

3200 Hz 


Theory 

Experiment 

Theory 

Experiment 

6.6 x 5.8 

0 db 

0 db 

0 db 

0 db 

4.4 x 2.9 

3.2 

2.8 

6.1 

-2.8 

2.9 x 2.9 

1.9 

2.0 

2.5 

-1.6 

2.9 x 1.5 

-0.6 

3.0 

3.1 

-3.3 

1.5 x 1.5 

-2.4 

3.4 

- 7.4 

-4.3 
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Table III. Relative SPL above Rectangular Areas of Plywood with 
Receiver Distance 2.7 m and Incidence Angle 70*. 


Scatterer 

Dimensions 

(m) 

£ * 

1600 Hz 

f - 

3200 Hz 

Theory 

Experiment 

Theory 

Experiment 

6.6 x 5.8 

0 db 

0 db 

0 db 

0 db 

4.4 x 2.9 

0 

-1.0 

0 

0 

2.9 x 2.9 

0.7 

-0.8 

0 

0 

1.5 x 1.5 

7.0 

0.8 

0.9 

3.4 


Table IV. 

Relative SPL above Rectangular Areas of Plywood with 
Receiver Distance 2.6 m and Incidence Angle 80°. 

Scatterer 

Dimensions 

(m) 

f = 

1600 Hz 

£ ■ 

3200 Hz 

Theory 

Experiment 

Theory 

Experiment 

6.6 x 5.8 

0 db 

0 db 

0 db 

0 db 

4.4 x 2.9 

2.3 

3.0 

6.5 

i 

O 

• 

tn 

2.9 x 1.5 

-1.6 

8.0 

3.1 

8.5 

1.5 x 1.5 

-3.4 

3.0 

-9.1 

6.0 




Table V. Relative SPL above Rectangular Areas of Plywood with 
Receiver Distance 2.4 m and Incidence Angle 80°. 


Scatterer 

Dimensions 

(m) 

f * 

1600 Hz 

f * 

3200 Hz 

Theory 

Experiment 

Theory 

Experiment 

4.4 x 2.9 

0 db 

0 db 

0 db 

0 db 

2.9 x 2.9 

0.6 

-0.1 

1.4 

-7.9 

2.9 x 1.5 

-0.2 

3.0 

1.0 

-10.2 

1.5 x 1.5 

-2.8 

0.6 

2.3 

2.0 


Table VI. Relative SPL above Rectangular Areas of Glass 

Fiber over Plywood with Receiver Distance 2.7 m 
and Incidence Angle 80°. 


Scatterer 

Dimensions 

(m) 

f • 

1600 Hz 

f - 

3200 Hz 

Theory 

Experiment 

Theory 

Experiment 

6.6 x 5.8 

0 db 

0 db 

0 db 

0 db 

4.4 x 2.9 

-0.3 

-2.6 

0.2 

1.0 

2.9 x 2.9 

-1.4 

-1.6 

0.1 

0 

1.5 x 1.5 

-2.6 

0.4 

0.1 

-3.4 
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Table VII. Relative SPL above Rectangular Areas of Glass Fiber 
over Plywood with Receiver Distance 2.7 m and 
Incidence Angle 70°. 


Scatterer 

Dimensions 

(m) 

f « 

1600 Hz 

f - 

3200 Hz 

Theory 

Experiment 

Theory 

Experiment 

6.6 x 5.8 

0 db 

0 db 

0 db 

0 db 

4.4 x 2.9 

-0.6 

-1.0 

0.1 

-3.0 

2.9 x 2.9 

0 

-1.2 

0.1 

-4.7 

1.5 x 1.5 

-1.0 

-2.8 

-0.4 

-5.4 


Table VIII. Relative SPL above Rectangular Areas of Glass Fiber 
over Plywood with Receiver Distance 2.6 m and 
Incidence Angle 80°. 


Scatterer 

Dimensions 

On) 

f * 

1600 Hz 

f - 

3200 Hz 

Theory 

Experiment 

Theory 

Experiment 

6.6 x 5.8 

0 db 

0 db 

0 db 

0 db 

4.4 x 2.9 

-1.1 

-10.0 

0.6 

-13.8 

2.9 x 2.9 

-1.9 

-2.5 

0.4 

-21.8 

1.5 x 1.5 

-0.8 

-7.5 

0.2 

-3.4 


Table IX. Relative SPL above Rectangular Areas of Glass Fiber 
over Plywood with Receiver Distance 2.6 m and 
Incidence Angle 70°. 


Scatterer 

Dimensions 

(m) 

f - 

1600 Hz 

f - 

3200 Hz 

Theory 

Experiment 

Theory 

Experiment 

6.6 x 5.8 

0 db 

0 db 

0 db 

0 db 

4.4 x 2.9 

-0.6 

5.5 

0.1 

-0.9 

2.9 x 2.9 

-0.9 

5.5 

0.2 

-0.8 

1.5 x 1.5 

-0.6 

6.0 

-0.4 

3.0 
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Appendix A 

Approximations for the Incident Pressure 
and the Great's Function 

It is desired to obtain an approximation for the integrand of 
equation (5) in which the source and receiver distances from the 
center of the scattering area are used as reference quantities, 
correction terms being incorporated in the exponent in the integrand 
to include diffraction effects and in the remaining factors in the 
integrand to indicate additional dependences on the size of 
scattering region. 

In order to accomplish this, it is expedient to expand the 
factors in equation (5) which involve the distances |r - r Q | and 
and | r^ - r^ | as power series in x Q and y Q , yielding (to second 
order) , 


r - r Q | » r(l - */r + V/r 2 ) 

(Al) 

r * rj' 1 : r _1 (l + */r - Tr 2 ) 

(A2) 


with 


- sine (x 0 cos* + y Q sin$) 


(A3a) 
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* • 

L ^v 


V(e ■ — (1 - sin 2 0 cos 2 *) + — (1 - sin 2 e sin 2 *) 
2 2 


1 .2 . „ 

T Vo sin 0 sm2 * 


(A3b) 


T(0,*) = — (1 


y 2 

3sin 2 0 cos 2 *) + — (1 - 3sin 2 0 sin 2 *) 


3 2 
-- x Q y Q sin e sin2* 


(A3c) 


Applying these approximations throughout equation (5) and factoring 
out the constants results in an integral of the form 


/■ 


-ikF(x ,y ) 

dS- e ° <?(x 0 ,y o ) 


(A4) 


in which the abbreviations are 


F - * + * s - (V/r ♦ V s /r s ) 


6 - 1 + 


(A5a) 


H - ?) * 0 a - * 2 ] 

[ (5aW i * -2*iBB s ® s ) 


(A5b) 
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and 


B 


1 + i(krB)* 1 ; E 


nB cose 
1 + nB cose 


E 


1 + 


B 

1 + n B cose 


(A5c) 


Upon collecting like powers of x Q , y Q and introducing the change of 
variables X * 2x^1, , Y = 2y 0 /W the expressions for F and G become: 


F = -(^X 2 + 6 2 Y 2 ) + (a 2 + yX)Y + a x X (A6) 

G = (1 + MX + NY + QX 2 + FY 2 + SXY) (n - n A ) (A7) 


with 


M = j sine cos<|»^ + (same) ; 

/W B \ 

N = l — — sine sin $J * (same) s 


I ± riii® 2 ...2, . I fl . vin 2„ cos 2jt 


Q ■ - sin 0 cos ♦ - - (1- 3s in 0 cos <j>)J j- 


if / \ \ 

+ Jsamej ♦ rF’\ E ^ E s^ 2+ J ^s^s] sine s * ne s cos * cos< *s 
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R = (1-EB) 2 sin 2 e sin 2 * - - (l-3sin 2 0 sin 2 *) 

Ir L 8 


']! 


( b 2 

/- - , 1 

— \ 

same} + 

(E+E-2 + - 

EE e BBj 

(s rr s 

\ s 2 

s sf 


sine si 


sine s sin* sin* s 


o (LW . 2 . „ 

S = < — sin e sm2* 




+ Isamel 


LW /- - 1 — \ 

+ 2rr~V E + E s ' 2 + 2 ^s^s / Slne sine s sin C* + * s ) 


and 


■ 7 k (w)( $i 


COS 


sme si „t * 


► . „ COS 1 

sine . * } 
s sin T s / 


6 


1 " 8rr 

2 s 


( l2 \\ ( 2 cos2 *\ 

\»2 )[ r s l 1 ' sin 6 sin 2 / 


* r 1 ' si " 2e s s& ♦* 


M 


(A9) 


Y * (r_ sin^e sin2* + r sin^e. sin2*_) 

8rr„ s s s 


These expressions are to be used in equations (6) -(8), 



Chapter 5 


PLANE WAVE DIFFRACTION BY A WEDGE 
WITH FINITE IMPEDANGE 


This chapter consists of a paper by Allan D. Pierce 
and W, James Hadden Jr. which appeared in the Journal of the 
Acoustical Society of America, volume 63, pages 17-27. 
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Plane wave diffraction by a wedge with finite impedance 

Allan D. Pierce and W. James Hadden. Jr. 

School of Mechanical Engineering, Georgia Institute of Technology. Atlanta, Georgia 30332 
(Received 8 September 1976; revised 30 July 1977) 

A theory is presented for the diffraction of acoustic waves by barriers with finite acoustical impedance, 
the shape of the barriers being such that, insofar as diffraction into the shadow zone is concerned, they 
may be idealized as semi-infinite wedges. The analytical development is based on the known exact 
solution for plane wave diffraction by a finite impedance wedge, versions of which have been previously 
given in the literature by Williams [Proc. R. Soc. London Ser. A 252, 376-393 (1959)], by Senior 
[Commun. Pure Appl. Math. 12, 337-372 (1959); Proc. R. Soc. London Ser. A: 213, 436-458 
(1952)], and by Maliuzhinets [Sov. Phys. Acoust. 1, 152-174, 240-248 (1955)]. This solution is 
described in detail and the asymptotic limit is derived in a form which demonstrates the satisfaction of the 
reciprocity principle. Practical implementation is discussed, both through numerical examples and through 
the presentation of graphs of quantities which will be helpful in barrier design. 

PACS numbers: 43.20.Fn, 43.20. Bi 


INTRODUCTION 

While the diffraction of sound around obstacles is a 
classic problem in wave theory, dating back to Poin- 
care 1 '* and Sommerfeld, 9 the design or assessment of 
proposed designs of barriers to reduce noise levels in 
areas adjacent to community noise sources is currently 
a topic of considerable interest in applied acoustics. 4 “ 7 
Ideally, such designs should be based on a comprehen- 
sive and accurate theory of sound diffraction around 
barriers. In practice, however, the inherent complexi- 
ties associated with the development of such a theory 
have necessitated the introduction of a variety of approx- 
imations and idealizations. Because of the remanent 
analytical difficulties, it is difficult to assess the ap- 
plicability of such approximations and idealizations to 
actual or proposed oarriers. In one of the most severe 
idealizations, the barrier is assumed to be perfectly 
rigid. Within the context of this idealization, it is prob- 
ably fair to state that the current status of the available 
theories and computation procedures is relatively satis- 
factory. 4 Diffraction around rigid barriers with planar 
surfaces can be considered using results derived from 
theories based on the ideal models of thin screens, 4,6 
wedges, 1 ' 9 and trapezoidal (three-sided) barriers. 7 ' 9 * 10 

The rigid-barrier theories, however, give informa- 
tion only on the effec s of barrier size, shape, and ge- 
ometry on diffraction; they give no insight into the ef- 
fect of the surface properties on sound levels in the 
shadow zone. Conceivably, the latter should be an im- 
portant consideration in barrier design. It is well 
known, for example, that the finite impedance of the 
ground may drastically alter the sound levels received 
near the ground from a source also located near the 
ground (i,e. , the so-called excess ground attenuation 
effect 11 ' 11 caused by the interference of direct and phase- 
shifted ground reflected waves). 

As regards available theories on the effect of surface 
Impedance on sound diffraction by barriers, the only one 
specifically devoted to acoustic diffraction of which we 
are aware is that of Jonasson 19 who gives an approximate 
theory of sound diffraction by a wedge of finite imped- 


ance. This theory, however, applies at best only to 
highly obtuse wedges, i.e. , where the exterior angle 0 
is only slightly greater than 180°. Moreover, it suffers 
from a lack of rigorous basis and is cumbersome to ap- 
ply: a crucial set of variables is presented only pictori- 
ally. Furthermore, a completely separate construction 
must be performed and several variables reinterpreted 
in order to show that the reciprocity principle 14 is sat- 
isfied (the point source solution should be invariant on 
interchange of source and receiver locations). It is ac- 
cordingly suspect, notwithstanding its good agreement 
with a limited amount of field data. 

There is, however, in the electromagnetic wave prop- 
agation literature , an exact solution for diffraction of 
plane waves by wedges of finite conductivity. Versions 
of this theory have been independently given by Wil- 
liams, 11-17 Senior, l9 ' 19 and by Maliuzhinets. 20,21 (Of the 
three, we have found Williams’s account 15 to be the most 
readable, although it suffers from a number of minor 
misprints and algebraic errors.) The purpose of the 
present paper is to extend and apply this theory to prob- 
lems of acoustic wave diffraction by wedges of finite im- 
pedance. 

I. STATEMENT OF PROBLEM AND SUMMARY OF 
RESULTS 

In this section we first describe the mathematical 
model on which our analysis is founded. Immediately 
following this statement of the problem, we present a 
concise summary of formulas for the estimation of the 
acoustic pressure diffracted around a wedge with finite 
acoustic impedance. This statement of results prior to 
their derivation is intended to facilitate the application 
of the results and to give an indication of the objective 
of the theoretical development in the following sections. 

A. The model 

We consider sources of such an extent and/or distance 
from the barrier's tip that the incident pressure waves 
may be approximated as plane waves. The geometrical 
arrangement is depleted in Fig. 1; the z axis of a cy- 
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PLANE 

WINES 



LISTENER 


FIG. 1, Diffraction of Incident plane wave by wedge of finite 
impedance. Listener coordinates are r. 6. a. The wave is 
incident from the 6 0 direction, wave-front normals make an 
angle > with the wedge edge U axis*. 


lindrical coordinate system is taken along the apex of 
the wedge; the surfaces of the wedge are the planes 
6 = 0 and 6-/3, where 0 is the exterior angle of the wedge 
(£>*)* We consider plane waves [with time dependence 
(e~* wt ) suppressed throughout the analysis] incident from 
the direction 0 O and at an angle > with respect to the * 
axis. On the surfaces 0 = 0 and 0 = £, the acoustic im- 
pedance is given in terms of a dimensionless quantity 
r) as 

Z *p 0 cv , ( 1 ) 

where p 0 r is the characteristic impedance of air. This 
description of the problem is amplified in Sec. II; it 
should suffice however, in the explanation of the nature 
of the results. 


AIL=-101og„{|l + |S^,e 0 )'(»jsin>)ll 1 } (5) 

in which one must use 

s/e , e 0 ) = 2 1 ♦ e 0 ) ♦ M v (e - e 0 ) |-‘ - qj, - 9) - qj,- e,) 

<e) 

with Q^(-0) obtained from_Fig. 2 or Fig. 3. Further 
discussion of the function QJi- 0) is presented in Appen- 
dix D. Several numerical examples, in which the com- 
putations may be performed using modern desk calcu- 
lators , are discussed in Sec . VII . The analytical steps 
which intervene between Pts. A and B of this section 
are discussed in the following sections. 

II. FORMAL SOLUTION FOR DIFFRACTION OF 
OBLIQUELY INCIDENT PLANE WAVES 

In the present section, the formal solution is summa- 
rized for the diffraction of obliquely incident plane waves 
by a wedge of finite acoustic impedance. This is es- 
sentially the same as those solutions given previously in 
the literature by Williams, 1 * by Senior, 1 * and by 
Maiiuzhinets, 10 although with considerable changes in 
nomenclature. Consistent with the discussion in Sec. 

I. A, the incident plane wave is taken in the form 

/>i, c = exp(- i>rsin> cos(6 - 0 o )|exp(i*z cos> ) . (7) 

Here 6 0 denotes the angular coordinate of the direction 
from which the incident wave is coming, > (taken be- 
tween 0 and i 77) represents the angle which incident 
wave- front normals make with the 1 axis; k is u?/c. One 
may note that the z -translational symmetry of the prob- 
lem implies that the resulting solution for the acoustic 
pressure should have the same z -dependent factor as in 
(7) above. The dependence on 6 and r is governed by 
the reduced wave equation 


B. Estimated insertion Ion 

For purposes of barrier selection or design, it is de- 
sirable to have an estimate of the effectiveness of a 
barrier in reducing sound levels at a given location. 

Ease of computation is certainly desirable. The model 
should be a reasonable idealization of practical cases. 
These conditions are fulfilled for the case of plane waves 
diffracted by a nearly rigid wedge with exterior angle 
£(>*) for larger observer distances r from the wedge 
tip, viz. , such that the condition krsiny » 1 holds (where 
w/c), and for angles 0 considerably less than 0 O - * 
(i.e. , listener well Inside the shadow zone). 

The quantity of interest is the insertion loss 

IL«20tog|o(|P MMr.l/li* witkMT.1) (2) 

which, in the case of a rigid barrier, is well described 
by the formula 9 

IL- 101og, 0 (*rsiny) - 201og l0 [.W;‘(» - *,) + M;‘(* ♦ * 0 )) 

(3) 

In which we have used 

'“‘■- •"fl ag ” <*> 

and v « ff/0. The principal result of this paper is that 
for a hard (but not rigid) wedge, there is an additional 
term in the insertion loss estimate, given by 


[ 8 T^ + ^ + 7 I ^ r + fr * 8inS ] /, = 0 ’ (8) 

Boundary conditions at the wedge faces are that the ratk> 
of pressure amplitude to inward normal fluid velocity 
component amplitude be p 0 rq, where tj, the specific (di- 
mensionless) normal incidence impedance, should have 
a real part greater than zero for an absorbing wedge. 
(Typically its imaginary part is positive, although not 
necessarily. ) Thus one has 

Bp /BO ± ( ikr/ii)p = 0 , at 0 * 0 and 0 = /? , (9) 

where the upper and lower signs correspond to 0 = 0 and 
0, respectively. 

An alternate parameter describing the wedge imped- 
ance which proves to be especially eonvenient is that 
of the (complex) angle or, defined such that 

C080T * (17 siny)" 1 (10) 

and such that - i *< < i it, a, > 0 given n„ > 0 and siny 

positive. The sign of a* is determined from sgn(oii) 
*sgn(T?|). For a rigid wedge, Fora 

perfectly soft wedge, tj- 0,0-1^. 

The solution for the boundary value problem as posed 
above may be taken in the form of a contour integral 





20 


A. D. Pitre* and W. J. Hadden: Plant wtvt diffraction try a wtdgt 


97 20 



FIG, 4. Integration contours in the complex { plane for evalua- 
tion of the acoustic field caused by a plane wave incident on a 
wedge of finite impedance. 


/> = exp(t*i cosy) ~ £ *xp(- 


ikr siny eos£) 

'/(t.M*.*)* , 


(ID 


where the contour C ( for the t integration may be taken 
(see Fig. 4) as C, +C„ + C„, where C, Is the path of 
steepest descents passing through the saddle point at 
t <0 of the exponential factor in the integrand, going 
from t = to £ = Similarly, t' n is the 

path of steepest descents going r rom £-aff + t*>to£ = j 
- 1 « through the saddle point at t = s. The contour ^ in 
encircles in the counterclockwise sense all poles of 
fit , 9 . <* 0 . “) which lie in the t plane between 0, and t . 
Since f (described below) is an odd (Unction of { , the in- 
tegral on contour c, vanishes identically, so only con- 
tours C|, and Tin are of interest. 


The function /(t, 9, ® 0 , a) is of a relatively complicated 
form and given by 


- S(t - 9 , *o. «) Wt - 9, 0 O , a) 


with 


... „ . (v/2)sln(vd a )*.(t, i * - a - p) 

Here we have abbreviated 

*„(a , h) - sinl (i •>)(« ♦ 6)) sln| ( J i >)(« - 6)1 

• Hc0S(t'ft)-C08(w>)j , 

V = ff/tf . 


( 12 ) 

(13) 


(H) 

(15) 


(Note that h is an even function of t.) The function S is 
defined by 


SU,9 9 ,a) = H^,a)/H^-9 n ,a) , (16) 

where the function ff.(t.a) is defined in terms of a (Unc- 
tion F/t) (here termed Wiliams' F function in recogni- 
tion of the fact that it is the same as used by Williams 1 *); 




F£t t£zJL±£!} Etii + t1f ) 

F t ({. + 2<? ♦ a) FJit + 2d - 1») 


(17) 


The analytic properties of the functions F££) and 
NJit ,a) are discussed extensively in Appendix A . For 
angles 0 of the form 


£=/>ir/2 q 


with /> an odd Integer and />, q relative primes, the func- 
tion FJX) is given by 

F^{)= n stn{| v|c*i»(4*»-l)-2^1}/ 

«fl ' 

•JI sin{i It — i w — + 1)1} . (18) 

Met 

Expressions for F^t) can be obtained for other values of 
£, but at the expense of considerably more computational 
effort. 15 


That Eq, (11) is indeed the appropriate solution can 
be ascertained by explicitly substituting it into Eqs, (8) 
and (9) followed by some integrations by parts. The 
fact that / is the sum of a function of £ - 8 and a function 
of t f 0 is sufficient to insure that the partial differential 
equation be satisfied. The boundary conditions are sat- 
isfied by virtue of the manner in which N $ (t,a) is de- 
fined in terms of Williams * F functions and of the tact 
that the h'& in Eq. (13) are periodic in £ with period 
20 . The explicit form of the function h was chosen in 
conformance with notions of radiation conditions; i.e, , 
that at large r the solution must consist of waves (other 
than the incident wave) which proceed outwards from 
the wedge and which do not grow exponentially with r. 
This requires in particular that / not have any poles 
between C\ and C„ for which the imaginary part of cost 
is positive. Since the function HJit - o) does not nec- 
essarily have this property, li(f: — or) was designed 
to have a zero which just canceled the "forbidden*’ pole 
of HJit - 8, or). Also, in order that the solution repro- 
ckice the assumed Incident wave, it was required that / 
have poles at £ = 8 - * 0 and at £ = 8 0 - B one of which is 
enclosed by C x „ when geometry indicates the incident 
wave is present. Finally, the function was required to 
have residues of appropriate values at these poles such 
that the C IU integration would give a term in the evalua- 
tion of (11) equal to (7) when geometry indicated the 
presence of the incident wave. It has been verified 11 
that this formulation is consistent with notions of reci- 
procity. 

The limiting cases of rigid and soft wedges may be 
obtained by examining the limiting rorms of the functions 
Ntf, a). S(C,0 o ,a), and In the limit of a rig- 

id wedge (o ’•lx), the limiting form of f(t , 0,fl Ol o), Eq. 
(12) is 

it)- 9,(t, 6-6*) + 9*9 t ) 

with 


(19) 


m 


* 


wm 
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QM , 0) a sin( i*) ' l cos(i/0) - cos( vt )] . 

Equation (19) corresponds exactly to the function used 
in Eq. (1) of a previous paper 9 by one of the authors. 

In the limit of an acoustically soft wedge (a -iao), the 
function / becomes 

/<t . Ms. * <?,(£ . » ♦ ®o> (20) 

which, again, is the correct limit. 

III. ASYMPTOTIC SOLUTION FOR DIFFRACTED 
WAVE 

In the shadow zone, the major contribution at large 
distances r from the edge comes from the C n portion of 
the contour integral in Eq. (11). The C UI integration 
simply gives the incident wave, a specularly reflected 
wave, and possibly a surface wave; the former two of 
which do not exist in the shadow zone, the third of which 
generally dies out exponentially with large r. The con- 
tribution /> Dlff# n from contour C lt at large r may be ob- 
tained by application of the saddle-point approximation 
taking into account the possible proximity of poles and 
zeros to the saddle point at £ = v. 

The poles of /(£ , 0,0 o ,o) are (i) those corresponding 
to the incident and specularly reflected waves and (ii) 
any pole of -0,a) or HJ-t -0,ar) which is not also 
a zero of *„(- £ - 0, »/2 - a - / 3) or *„(£ - 0, ff/2 - a - 0), 
respectively. [See Eqs. (12) and (13), ] The first cate- 
gory of pole is manifested by the factor *„(0 O ,£) in the 
denominator of the definition (13). The second category 
of poles may be determined with reference to Eqs. (A6); 
the only ones which could conceivably be close to the 
£ = 7? saddle point are where £ - 0 or - £ - 0 equals f tt 
t a or - l it ± a f 0, respectively, or thus where £ = \ it 
* a + 0 , £ = §ff±a + 0~0. For given 0 and a , at most 
one of these poles will be near the saddle point. Let 
us assume that the relevant pole is at £ = 7 r + ; we then 

set 


= . 


( 21 ) 


Consequently, if one sets 


/(£,M 




( 22 ) 


where, upon rearranging the product of **(0o»t -#) and 

*„(0o»t + 0) into the f actorB 


D 4 «cos(W:)-cosK0-0o) , 

Dj = cos(y£)-cosi>(0 + 0 o ) , 


(23a) 

(23b) 


The function *(£, 0, 0<>, or) so defined will have no poles 
in the vicinity of £ = v, 

T>e analysis then proceeds, as described in Appendix 
C, by replacing 4>, Z>j, £>,, Z), by power series expansions 
to first order in (£ - ir) and integrating the resulting form 
of Eq. (11) along the line of steepest descents through 
the saddle point at £ * s. Thus the integral on contour 
Cj, in Eq. (11) can be expressed in terms of standard 
functions occurring in diffraction problems as 


n =exp(f*(zcosy + rsiny)](e' #/4 //50 


y [ g 1 "a „< rMl **) + c‘ -'A oirAC') + 

c-u^rpjl 


(24) 

where 


^(♦» G(ff,0,0 o ,a) 

(25) 

A/j->=Af„(0-e o ): Afi‘* = Af„(e + e 0 ) , 
with M'„** given by Eq. (4). and 

(26) 

G(£,e,e 0 ,a)«+(£,M 0 ,«) (rslnrr) 1 . 

(27) 


The other two coefficients of the function A D (X) in Eq. 
(24) are obtained by cyclic permutations of the quantities 
A AfJ‘\ and P* . In the argument of A D (X), we have 
used 

r z[(kr siny) . (28) 

The function Ap(X) which appears in Eq. (24) is the dif- 
fraction integral defined in the previous paper 9 by 

W r e-*ds 

Ad(x) u J . v29 ' 

which, when X is real, can be expressed as 

= sign(X) ( /( i X| ) - | xj )] , (30) 

where /(X) and g(X) are the auxiliary Fresnel functions 
tabulated on pages 323-324 of the NBS Handbook of 
Mathematical Functions ** If Xis not real, as would be 
the case for X- TP a , the above would be inapplicable, 
but, instead, one could write 

A 0 (X) = (<?*'" 4 Nl) u{e l1 "'(ii/z) xlt x] (31a) 

or 

ylpU) = - {*'*'* /&) V 2) ui X* I* (31b) 

which would hold for Im(e ,,/4 AO positive or negative, re- 
spectively. The function w(z) is related to the error 
function of complex argument and is tabulated on pages 
325-328 of the NBS Hattdbook 23 for complex values of z. 


The calculation of the function C(£, 0, 0 O , <*), which is 
quite tedious, is sketched in Appendix C. The results 
which are relevant to Eq. (25) are 

C (M, 0 O . Of) = P. ^0 O , a) U{9 f ^)/?(0, 0 O , nr), (32) 

where 


l t0,ar) 


($)atn(i/fl) 

tf,(-e,or)*.(^,il»-or-0) 


(33) 


and 


^0 f 0 O , <*) = Afv(0+ 0 O ) + AU0 - 0 O ) 

cos(2i/Qf) — cos(uff) 

’ ysim'ii 


(34) 


In the complete asymptotic limits where P a , A/i°, 

M[ m) are all finite and T is large, A D (X) can be replaced 
by (*X)* 1 and a considerable simplification results in 
Eq. (24). Specifically, one finds 

/Wjr.11 =exp[»(«co«r+ rainy)] 

y («?' ,rt /i/5) (try'G(0, 9 0 , a), (35) 


where 


* 


* 


A 


1m .1 i 
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« 


G(0, 0 Ol a) * C (*, 0, 0 # , o)/P. Af i*‘ Af J*» . (38) 

Note that, with G(i, 0, 0 O , a) given by Eq. (32), the (ac- 
tor P, cancels out. 


Making similar expansions In Eqs. (33) and (34) the 
function <5(t, 9, 9 t , a) in Eq. (32) may be written as 

G(w, 9, S 0 , o) * P.l Af * * * M, *| 1 1 ♦ $,{9, 0 O ) (i f - or) i (43) 


The symmetry of Eq. (36) with respect to 9 and 0, Is 
obvious from Af„(0 - 0 O ) « M,(0 O - 0). Thus the reciprocity 
requirement Is definitely satisfied. In the limit of a rigid 
wedge (a =$») one has 0(0, i*) « - 1 by virtue of Eqs. 

(14), (17), and (A3), and 0(0, 0 O , l w) - Af„(0 ♦ 0 O ) + Af„(0 - 0 O ). 
Consequently, one has 


G(0, 0„ i») 


1 1 

Af„(0 + 0 O ) + Af„ (0 - 0 O ) 


(37) 


and the result for asymptotic diffraction by a rigid wedge 
is recovered. 


and consequently, from Eq. (36), 


C(0,0o,o) 
in which 


* iOfH) ' * - ol 1 


(44) 


S,(0, 0o> * 2lA»,(0 ♦ *,) + Af „(0 - t* 0 ) I* 1 - Q,(- 0) - ^,(- 0 O ) . 

(45) 

For observation angles other than 0 *0, the diffracted 
pressure field may be approximated— by combining Eqs. 
(24), (25), (40), (43), and (45) -as 


Similarly, in the limit of an acoustically soft wedge, 
i. e. , o - 1 * *>, 0(0, e 0 * <») approaches exp(- «2vor)/ 2v simnr. 

i t - o - £) approaches - i explMw^ - a) I, 

W 8 (- 0. a) approaches exp| - (i*tt)], sot/(0 o , OfH/(0 o , a) 

x 0(0.^o» a ) approaches 2sinytfsind p divided by I'Simnr 
or just Af„{0 + 0*) - Af„(0 - 0 O ). Consequently, one has 


G(0, 0 o ,i*>) 


1 

A/„(0 - e 0 ) 


1 

A, M (0 + 0 O ) 


which, again, is the correct limit. 


(38) 


In addition to the above-mentioned contributions from 
the integration contour C n , a complete description of 
the pressure field in the shadow zone should include the 
possibility of a contribution from a surface wave which 
is refracted from the shadowed face of the wedge. Such 
a contribution would arise from a pole enclosed by the 
contour C m . For tK ir and 0 O > r + 0, the only pole of 
fitf 0* 0q» a) that could conceivably lie within C, u is at 
£*!*-<»♦ 0. This pole will lie within the contour only 
if 


i ff - < sin* l (tanhof t ) (30) 

and since in this case the imaginary part of cos(| * - a 
+ 0) is negative, by (11) the C ul contribution from this 
pole decays exponentially with distance from the sur- 
face. Furthermore, reference to Eq. (10) indicates 
that the inequality (30) is not likely to be satisfied in 
situations of physical Interest. For these reasons we 
omit an explicit description of the surface wave contri- 
bution (which is included in Sec. IV of Ref. 22). 


IV. NEARLY RIGID BARRIERS 

If the barrier is nearly rigid, for n^O a is close to 
i ir and one can take the solution of Eq. (10) as 

Jtr - a * fa siny)* 1 (40) 

Thus it would seem appropriate to expand //,(£, or) in a 
power series in J it - o, keeping up to first-order terms 
la i ir — or. In this event om has, from Eq. (17) 


or) * - UuiUpcIII ♦V*(t)(J*-of)|, 


where 




® + d ♦ 1 * )/**(£ ♦ 2*1 + | tr)l 

*1 ? -"f* ^(7T27rr| 7) J • 


(41) 

(42) 


Pdxut, u * explf M* cosy + rsin> ) I (c ** f 4 A 7 "?) 

* p. Uf l * 1 ♦ Af J* M 1 r '* U „(r \i ? >) 

V ,a U p (rRj| , (46) 

where 


jrl*) 


1 

[ A/J‘ ’ - A/ j* 1 1 ( P„ - A/ j* ’] 


(47) 


and the other coefficients in Eq. (46) are obtained by 
cyclic permutation of A/*\ A fj"\ and Similarly, in 
the complete asymptotic limit, one has 

/>Ditfr, !Vo*i ™ «xplf*U cos > Sin) )| (o u/ Wll) (nr)- 1 

♦ [Af t } Y x } ll ♦ 5,(0, 0 O >/*) siny] , 

(48) 

From this last expression one may obtain a correction 
to the insertion loss for a hard barrier {vis -a -vis a rigid 
barrier) given by 

- 10 log l0 | 1 4 S,(0, 0 o )/rj sin> | 1 dB (49) 

which is just Eq. (5). 


V. PRACTICAL APPLICATIONS 

The formulas presented in the preceding sections 
probably appear more formidable than is actually the 
case. The first important point is that these results 
can be used most fruitfully in the computation of barrier 
insertion loss, as set forth in Eq. (2) 

lL fcr . »201og,o(|/>„ 

Mir. I/I p tar* i ) • 

If, as is generally true, the surface wave contribution 
may be neglected, Eq. (5) provides an estimate of the 
change in insertion loss, with respect to a rigid bar- 
rier's effect, of a barrier with finite impedance. 

The second important consideration concerns the func- 
tion V»(- 0) which appears in this correction term, it is 
shown in Appendix D that for angles ^ of the form /nr/ 

2?, Q $ (- 0) has a form which is amenable to numerical 
computations; in some cases computations are not so 
taxing as to require a large computer. In practice, it 
should be possible to obtain useful estimates of the in- 
sertion loss (or sound pressure distribution) using a 
value of ^ of the above form. One then has [see Eq. 

(I>4)| 
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CONFIGURATION 12 3 4 


AlL(dB) 2 3 17 M 05 

FIG. 5. The finite -impedance correction to the ineertion loss 
for a wedge with interior angle 10 . The surface admittance 
was taken as (1 V) =0. 1 - <0. 05. The source -receiver orienta- 
tions are identified by the configuration numbers: in configura- 
tion 1 the incidence direction is at 30" from the adjacent wedge 
face while the receiver is at 45 6 from its adjacent wedge face. 


_ o-i)/* j 

Q a ( 0)- i/ sin(x/ir) s ^| v ^V2«»)Tsin{p[d -(2«- 1)*]} 

V' gin(6) + 2mg) + s in|0 + (2m + 1)0] . . 

+ sin(0 + 2w0) sinfeT + (2w + 1)0T 


As an alternate aid to applications of these results, 
we have computed Q a (- 0) for a number of values of 0 
and 0. These are presented in Fig* (2). In addition, 
these curves are plotted again, with 0 appearing as the 
independent variable, 0 as a parameter, in Fig. (3). 
Thus one has the option of using one of the “special” 
values of 0 to approximate a desired wedge or of using 
values for the desired wedge angle for a selection set 
of angles 0. 

As an illustration of the use of these results we have 
calculated the correction to the insertion loss via Eq. 

(5) for perpendicular incidence (ysjir) on wedges with 
exterior angles 0=350° and 240° for a surface admit- 
tance if 1 =0. 1 - i 0.05, which is representative of turf 
at 1000 Hz . n Values for the function Q were obtained 
from Figs. 2 and 3. Insertion-loss corrections are 
presented in Figs. 5 and 6 for several incidence and 
observation angles. 

In the case of the acute-angled wedge, the surface im- 
pedance has a small effect— less than 3 dB. For the 
obtuse-angled wedge, the effect from considering the 
finite acoustic impedance is on the order of 6 dB when 
both source and receiver are at fairly small angles with 
respect to their adjacent sides of the wedge. Thus con- 
sideration of the finite impedance seems to be of con- 
siderably greater significance for obtuse wedges than 
for acute ones, especially since in many practical real- 
izations of the obtuse wedge model the sources and re- 
ceivers are close to the surface. 


For angles 0 of the form 0-W2?, k an integer less 
than or equal to p, there is a singlular term in each sum 
in Eq. (50). A straightforward expansion of the two 
terms reveals that they combine so that Q a (- 0) is in 
fact regular. Difficulties in numerical computations 
may be avoided by avoiding such angles. 

A more detailed investigation of the cases in which 
the receiver angle 0 is very small or the source angle 
0 O approaches 0 reveals that Q a (- 0) becomes 

Qj(- 5) *<?*(- 0 + 5) » (sinfl)" 1 , 5<^1 . (51) 

Since § f (- 0)/(qsiny) serves as a first-order correction 
term to the rigid-wedge limit for J/ a (- 0) [see Eq. (41)], 
the behavior of Q $ exhibited in Eq. (51) indicates that 
the approximation in Eq. (41) is not useful for situations 
in which the incidence or receiver directions are at 
small angles 5 with respect to a wedge face. This be- 
havior is a manifestation of the familiar phenomenon 
of the vanishing effective surface impedance for plane 
waves at grazing incidence . u This is borne out by in- 
spection of Eq. (A5): Substituting in the appropriate 
values for £ and or, one has 

/?[$, i v - 0? siny)" 1 ] [0 - 0. 1 * - (q siny)’ 1 ] 

. ( tan(j i/5) - tan(v/2t) siny) j 
| tanti + tan(v/2f? siny ) ) 

for the plane wave reflection coefficient at each face. 
For a given value of 17, the reflection coefficient ap- 
proaches - 1 as 5 goes to zero. These considerations 
indicate that useful estimates of the insertion loss cor- 
rection can be obtained for 5r? siny » 1. 


VI. CONCLUSION 

A theoretical analysis has been presented for the dif- 
fraction of plane waves by a wedge of arbitrary surface 
impedance. Particular attention has been given to the 
pressure field in the shadow zone for large distances 
from this tip of the wedge. The results presented here 
make use of simplifications that result for a large num- 
ber of special wedge angles. In a detailed discussion of 


RECEIVER 
^ 3.4 



SOURCE 


CONFIGURATION 12 3 4 


AIL (dB) SI 4 6 IS 6S 

FIG. 6. The finite- impedance correction to the insertion loss 
for a wedge with interior angle 120’ and surface admittance 
(1 n) ~0. 1 * /0. 05. Source and receiver orientations are 
labelled as in Fig. 5. 
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the nearly rigid wedge a correction to the insertion loss 
of a rigid wedge has been obtained* Numerical computa- 
tions indicate a significant effect of finite surface im- 
pedance for obtuse wedges with source and/or receiver 
at a fairly small angle with respect to the plane of the 
wedge face. In addition, the solution presented here 
should provide a good point of departure for the analysis 
of diffraction by spherical waves or by broad barriers. 
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APPENDIX A: DISCUSSION OF THE F AND H 
FUNCTIONS 

The F function F a (£ ) appearing in Eq. (17) *s defined 
such that it satisfied the two recurrence relations 


F 4 (£+20)=F 4 (£)tan[&-l v)| . 
F 4 (£ * 20) = F 4 (- £) , 


(Ala) 

(Alb) 


and such that it is analytic and has no poles or zeros in 
the strip 0< £*<20, the function asymptotically approach 
ing zero as £j-± 00 . Explicit expressions for this func- 
tion for some particular values of 0 are given in Appen- 
dix B. 

Since F fl <£) has no zeros or poles in the strip (0,20) 
it follows from Eqs. (Al) that its zeros must be at 


£ 30* (| 7T + 0 + /O 

while its poles are at 
£ =0±(iir + 0 + /fl , 


(A 2a) 


(A 2b) 


where K * 2ttir ♦ 2m 0 ^ 0, n and m being arbitrary non- 
negative integers. (Of course, there is the possibility 
that a pole location may coincide with a zero location, 
in which case the function will have neither a pole nor 
a zero at the point in question. ) Examination of the 
locations of the poles and zeros of the function {sin[(|y) 
x (£ + \ rr)] f 9 (C + ff)}"\ reveals that these are identical to 
those given above, so one may infer that F*(£) has the 
property 


F 4 (£)F 4 (C + »)- 8lnI(i|/)( * t+iir)I , 


(A3) 


where A $ is some number Independent of £ (the precise 
value of which is immaterial). It then follows from this 
and Eq. (A2b) that the asymptotic values of F f (£) should 
be 

*V£)-A,e ,wm , £,-- <A4a) 

-M**- 1 *" 4 *. £,--«. <A4b) 

Analogous relations may be deduced for the Ug(C t or) 
function starting from the definition of Eq. (17). It sat- 

l 


isfies recurrence relations of the form 

ff«(£. a) m _ f sin(y£)-cos(yq) 1 
/f*(- £> a) * L9ln(v£) + cos(i/a) J 

tanj(j v)(£- 2* + or)l 
"" tan|(J v ) (£ + 3 w~- or) 1 

£ — 0, or) n /w _ \ 


(A 5) 


where R(£, a) may be interpreted as a plane wave re- 
flection coefficient. The zeros of ff d (£, or) are at 


£ =- 0* a± (3 jr + 0 + A) 


and 


£ = ± (it — or) ± (g tr + 0 + K) 


(A 6a) 


(A6b) 


any sign combination being a possibility, while its poles 
are at 


£ = - 0± a ± (2 ff + 0 + K) 


and 


£ =± (* - a)± ( i it + 0 + tf) 


(A7a) 


(A 7b) 


for, again, K = 2rtir + 2m0, w^O, m^O. Also, it follows 
from Eq. (A3) that 

M A (t f a)H 9 (C + it % a) 

* si»Mii£ ± i> 1 ~ aJ h (as) 

cosf^lstt- Jr) + a ]}cosU tH|(C+ 3u)-aj} 

APPENDIX B: THE FUNCTION F^f) FOR 
PARTICULAR VALUES OF THE WEDGE 
ANGLE 0 

The function F,(ir) defined in Eq. (32) of Williams' 
p:»per is fails to exhibit the proper asymptotic behavior 
in the limit as <*-»>«, e.g. , it does not obey our Eq. 
(A4). Accordingly, we describe here the construction 
of our function F,(C) for 0-pt/2q, with P an odd integer, 
q an integer. We shall Include two examples for wedge 
angles of particular physical interest. 

We begin by recalling from Eq. (A2a) that the function 
F t (C) has zeros at the values * 

+ + (Bl) 

which, for ^particular values of fi under consideration 
here, may be written as . 

(2 q/n) (£ - pv/2q) = ± (/> + 3$ + 4 nq + 2 tup) ■ * F t (n, m). (B2) 

Similarly, F*(£) has poles at (see Eq. (A2b) 

C-P*±(iv+P+K) (B3) 

which In turn may be written as 

(2 q/ f ) (£ - pi/2q) *±(p+q+ 4 nq + 2 mp ) ■ t F # («, m). (B4) 

We may use Eqs. (B2) and (B4) to represent the function 
F 4 (0 schematically by 


F|(£)» fl [*-F,(n, im)][x+F.(», m)] / fl [*-F 4 (n, m)][x+F,(», »>)) 

m,Rwb / 


(B5) 


3SX.'-ilia 
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There are several relationships between the pole and zero locations, F t (n, m) and F,(n, m), which are important in 
the development of F t (C). They are 

(B6) 


M-q)*F p (n,m), 

F t (n , m) = JJ(-» + |(p- 1), -m-l-g), and 
Fp(n + p. m - 2q ) = F 4 (m, m) . 

Substituting Eq. (B6) into Eq. (B5), with appropriate changes of limits, yields 

F,(t) = ii m*- F.(», w))[jf+F > (n,m))/l [ ([ l*-F,(n, m)l(*+F 4 (»,m)]. 

»a| f ItiO MeO 

After canceling out the common factors, we have 

-i • , + «-i 

F.(C) = II II l*- J>(n, m)][x+F # (»», »»))/ J J | J l*-F,(», m)|[v+F,(n, m)\. 

We next apply Gq. (B7) to the [* + />(», m)l factors in Gq. (BIO) and rearrange the factors to obtain 

-l • M -**«♦!> / * !/*<#-!* -(#♦!) 

F,( 0* U II U-F # (n,w)| U II l*-F.(«,»»)]/ II H[*-F,(n,m)] fl H [v- F*(n, »«)] (Bll) 

ns-DSiM) ■»-« / it.0 «-0 R— - <n4( 

which becomes, after use of Gq. (B8) and further adjustment of the product limits, 

-i • -i -l j • «-i *-i 

= II M II U-F,(w,m))/II Ht*" F,<«, m)] [] II l*-F # (w, m)]. (B12) 

He-l /2(*-l)"“« *••1/ *<#•!) mm-~ I ft«0 *-0 ft* -• WeO 

Note that the numerator of Eq. (B12) may be consolidated to read 


(B7) 

(B8) 

(B9) 

(BIO) 


•l J* j 4 f 1 

Num - n ni« - F p (n, W)1 / n m* - Fp(n, W)1 . 

n.-l/to-l) «.-* I «-W *<»-!> »*0 


(B13) 


Similar treatment of the denominator of Eq. (B12) yields a form similar to Eq. (B13), but with the limits, I n I =», 
0* m £ q - 1. Thus we have the final schematic representation for F*(t) 


F a (D = ft ft I*- F>,m)]/n fit x-F^n.m)] . 

Be* 1 /{(#•!) f MiO *••• 


A specific functional form for F 0 (C) which satisfies all the requirements of zero and pole locations is 

-1 , «•! V 


F.(t)= ii sin(v/2/>)[x-/>- (4n+l)$|/ ft 4^ 1 x ~P~^ m *^Ph 

n*.l/t(^4) I BteO 


(B14) 


(B15) 


where we have used the definition of F 4 (m, m), Eq. (A4). Finally, making the substitution, x = (2q/n)(C-pv/Zq), and 
simplifying, we obtain 

»/»(»-!> 

(B16) 


t / ) / H 

F 4 (C)= JJ slnl^/^C-ir + ^g^pHdw-Dl/IJ »i n lit-I»-(/>v / 2 ?)(m + l)|. 

I 


We conclude by quoting two examples for particular 
wedge angles 0 which are of physical interest. First 
for the right-angle wedge, 0 = Eq. (B16) yields: 


APPENDIX C. ASYMPTOTIC APPROXIMATION 

OF /(r.Mo.a) 


f t rl 

F»»/,(t) - gin(If)+cos(it) ' 


Secondly, for the oblique wedge, 0 = 5*/4, we find 


e, /Vi vT[cos(i:)-cos(fs) 


)] 


Other wedge angles may be treated with greater effort. 
It is noteworthy, however, that the expression we have 
obtained for F*(t) is considerably simpler than that ob- 
tained by Williams. It may be verified readily that Gq. 
(B16) exhibits the correct asymptotic limits prescribed 
by Gq. (A4). 


In obtaining an approximation for the function /(£ , 0, 
0 Of a) in the vicinity of the saddle point at £ = ir, the 
( B1? ) quantities D lf D Zf D % and *(t,0,0 o , a) in Gqs. (23)-(25) 
are expanded in powers of (t - t) to yield 

/(£,e,0„ o) = [r. - (f - ,)) lM: i » - ft - *)) - (t - *)] ’ 

(B18) (Cl) 

where we have used A/J*' = M„(6 + 0,) and Af‘*’ = Af„(0 - ?,), 
with Af„ given in Eq. (6), while 




(C2) 


Is assumed to b« expanded In a power series in £ - * up 
to first order in ( £ - »). 
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As regards the actual integration, one replaces cost 
by - 1 + (|)(f - nf In the exponential factor of Eq. (11) 
and integrates along the line of steepest descent of the 
approximate integrand, i. e. , along a line going oblique- 
ly downwards making an angle of 45° with the real axis 
and passing through the saddle point at *. The integra- 
tion variable is charged to s, given by 

(£-»)* (2/*rsiny) ul se m>,/ * 

the s integration now going from -*> to «. The approx- 
imate factor for / is also replaced by means of the al- 
gebraic identity 


It serves as a corection term in the function /f,(f ) for 
nearly rigid wedges. Since the function F t ([) takes on a 
simple form for angles $ of the form pw/2q, it is rea- 
sonable to expect that Q t (t) might also be cast in a rela- 
tively simple form in the same instances. 

One begins by noting that from Eq. (B16) one has 
!.'**«» 

In F,(f) = £ Msln{ii'(t -20+ii(4n- 1)]}) 

ml 

- £ ln(ain{| (f -i»- 20(m + 1)]}) (D2) 


tv+xs 1 (w + x a \ . . 

(a-s)(b -s){c-$) ~ {6-a)(c-a)\ a-s “*) 

* two additional terms obtained by 
cyclic permutation of a, b, and e. 

This yields the asymptotic approximation to the dif- 
fracted pressure field, Eq. (24). 

It is now necessary to obtain a suitable expression for 
6(f ,0, 9q, a). By comparing Eqs. (12) and (22) and their 
associated definitions, one may obtain an explicit rep- 
resentation for <t> 




X [*„(9o. { ♦ 0)0i - *>(00. £ - 0)0,] , (C4) 

where 


0ij E +>(*£-0.s»-a-/3)/f,(±£-0,») • (C5) 

Now, by using Eqs. (A8) and (14) in Eq. (C5) and per- 
forming several trigonometric manipulations, one may 
obtain 

\ 4 fl,(±£ -0*», or) *„(t ?0-*,$v- a- 0)- (C8, 

Then after substituting Eqs. (C6) into Eq. (C4), the re- 
sult into Eq. (C2), and expanding and recombining the 
terms 

*>(0,. f + 0)**,(£ - 0 - i », ar) - *,(0„ a - 0)*„(£ + 0 - i i , o) 
and using the definition 

’ a) */f,(-0, at*>(0,i»- or-d) (C7) 

one has finally 

C(r,0,0„ or)- J>.C/(0, a)U(6 t , a)D(8,9 t ,a) (Cfl) 

with 

w. ... .) ■ m.(» . «.i . m.(» - ,.i . 

(C9) 

which completes the outline of the analysis leading to 
Eqs. (32)-(34). 


APPENOIX D. THE FUNCTION Q 0 iS) FOR 
PARTICULAR ANGLES 0 


The function Q t ({) is defined in Eq. (42) as 

WC) H ta ~ % _V»)V( £ * 20 - T *) J ’ 


(Dl) 


and thus 

t 

-jr lnf,(£)= | v £ cot{|vlt-20+|ir(4»-l)]} 

n-1 

1 rl 

- cottilf-lw-Zutm+D]). (D3) 

Then upon substitution of Eq. (03) into Eq. (Dl), con- 
solidation of arguments of the several cotangents, and 
use of the identity 


COt0j ± cot(0 2 ) * 


sin(0 2 *0 1 ) 
sin0, sin0 2 


followed by the use of trigonometric angle-addition re- 
lationships, one may obtain 

§,(£) — vslnM g sinlv(r + rnvTTsinMt" - 1)]} 

_ T s >n(t - 2wfl) * sin(£ - 0(2 m .1)] . . 

” Zi sin(£ - 2 m 0 ) sinlf - 0 ( 2 m +1)] ' ' 

It should be pointed out that for £ = - k/2q, with k a posi 
tive integer less than p, there is one singular term in 
each sum in Eq. (04). It can be shown, however, that 
the two singular terms combine in such a way that Q(£) 
is continuous at the apparent singularity. For f = 0, 
there is a true singularity in Q t . In this case one may 
see from Eq. (41) that this singularity is cancelled by 
the tan(iv£) factor in H t to give 


for a' *i*-(T)Siny)~ : , which indicates the manner in 
which the rigid wedge limit of HjiO, a) is approached as 
the impedance becomes infinite. Similar behavior may 
be noted for t » - 0- 
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Chapter 6 


EFFECTS OF AMBIENT FLOW AND 
DISTRIBUTED SOURCES 


Diffraction in the Pressure of Ambient Flow 

The model sketched in Fig. 1 may be used to assess the 
effects of ambient flow on barrier diffraction. The barrier 
is taken as a thin screen which occupies the region x < 0 of 
the y = 0 plane; the edge of the screen lies along the z axis. 
The source is taken as being localized as a point *g,yg,Zg 
where yg c 0, the listener is at (x^,y^,z L ). A uniform ambient 
flow of velocity U Q is in the +x direction, tangential to the 
screen and having the same velocity on both sides of the 
screen. Since the screen is idealized as being arbitrarily 
thin, there is no discontinuity in U Q at the edge. 

The solution for plane wave diffraction in terms of such 

1 2 

a model has previously been given by Candel. * Here, a 
slightly different approach is used for the case when the inci- 
dent wave ensues from a point source. 

If one limits one's consideration to a single frequency 
component and uses the device of taking e"* wt to describe the 
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Fig. 1. Geometry used in the study of the effects of 
ambient flow on sound diffraction. A point 
source is immersed in a uniform steady flow 
near a thin screen. 


time dependence, then the complex amplitudes of the acoustic 
field variables satisfy the equations 

[ - ia> + U Q 3/3x]p/c 2 + P Q V*u = 0 (la) 

P 0 [-iw + U Q 3/3x]u + vp ** 0 (lb) 

everywhere except in the immediate vicinity of the source. 

From these one may derive the generalization of the scalar 
Helmholtz equation which takes ambient flow into account, i.e.. 


£ x {p} = 0 (2) 

£ x = v 2 - c' 2 [-iu> + U 0 3/axl 2 (3) 


There is a transformation^ which, for U Q /c<l, will reduce 
equation (2) to one resembling the scalar Helmholtz equation 
without ambient flow, i.e., 


> iP . e _ i (M«/ c) x/B 2 £ J p e i(Mw/c)x/3 2 | 


(4) 


where 


M * U 0 /c 


( 5 ) 


is the Mach number, and 
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8 - [1 - M 2 ] 1/2 (6) 

£ 2 = a 2 /a(x/e) 2 + a 2 /ay 2 + a 2 /az 2 + (u>/e) 2 /c 2 (7) 


Consequently, one may conclude that any solution of Eq. 
(2) corresponding to a given angular frequency u may be taken 
as 


p(x,y,z,u>) 


e -i(Mw/c)x/8 2 


p(x,y,z,u) 


( 8 ) 


where 


x * x/8 (9a) 

y * y (9b) 

z ■ z (9c) 

u = u/e (9d) 

and 

[a 2 /ax 2 + a 2 /ay 2 + a 2 /az 2 + C^/c) 2 ]p = o (10) 


The latter is the scalar Helmholtz equation corresponding to 
no ambient flow. 

For the screen diffraction problem, one requires that 
Uy ■ 0 at y * 0 for x < 0. Consequently, from Eq. (lb), 3p/ay 
should also be zero for the same circumstances. However, one 
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sees from Eq. (8) that this requires 

3p/3y =0 y = 0 , x < 0 (11) 


which is the same boundary condition as one would have in the 
absence of ambient flow. Also, if the source is at Xg,yg,Zg, 
then p should correspond to a field generated by a source at 
x s ,y s ,£ s . Consequently, one concludes that the solution of 

£jP = -4irS(x,y ,z) (12) 


subject to the boundary condition mentioned above may be taken 
as 


P 



i(Mw/c)(x -x)/B 2 

e s ( x 0 ,y 0 ,z 0 )G(x/e,y,z| x Q /e ,y Q , z Q | w/e) dx o dy 0 dz Q 


(13) 


where G(x,y,z|x 0 ,y 0 ,z 0 |u) is the Green's function for the 
scalar Helmholtz equation in the absence of ambient flow, 
satisfying 

[V 2 + (w/c) 2 ]G(x,y,z|x 0 ,y 0 ,z 0 |a») - -4*6(x - x 0 ) (14) 

In the case S(x,y,z) is taken as S Q 6(x - Xg) , the above reduces 
to 
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p 


i(Mw/c)(x s - x)/e 2 

e G(x/8,y,z|x s /&,y s ,z s |w/e) 


(15) 


which may be considered as the Green's function for the 
problem of diffraction of waves by a thin screen in the 
presence of ambient flow. 


Effect of Ambient Flow on Insertion Loss 

The Green's function without ambient flow included is 
described in some length in Chapter 1 of the present report 
and the fact that it is amenable to numerical computation 
implies that the problem discussed above is also. Here, for 
simplicity, we limit our discussion to the circumstances 
described by Fig. 2. The Green's function in the absence of 
ambient flow for such circumstances is given by the Fresnel 
number approximation 


e ikL e^/4 

L ST 


f([2N) 1/2 ) - ig([2N] 1/2 ) 


S. 

L 


ikL 


(1/2) - N 1/2 e‘ lir/4 


N << 1 


ikL iir/4 
e e 


2ttN 


1/2 


N >> 1 


(16a) 

(16b) 

(16c) 


where N is the Fresnel number given by 


Ill 



Fig. 2. Sketch of limiting case 


discussed in the text. 
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N 


L - R 
(x/2) 


(17) 


where x = 2n/k and 


[(z L - z s ) 2 ♦ 

( r L * r s ) 2 ) 1/2 
# 

(18a) 

[(x L - Xjp 2 + 

(y L - y s ) 2 * (* L - V 2 I 1/2 

(18b) 

C 4 ♦ ^ 1/2 *> 

r s ■ < x I * yl) 172 

(18c) 


The functions f and g are the auxiliary Fresnel functions 

4 

tabulated in the NBS Handbook of Mathematical Functions. 

For most purposes, the asymptotic limit (16c) may be con- 

1/2 

sidered as realized when (2N) > 2 or N > 2. 

The insertion loss due to the barrier is defined as the 
loss in decibels of the sound pressure level at the listener 
location due to the presence of the barrier and, for waves 
from a point source, is accordingly 

IL - 10 log 10 |G NB /G B | 2 (19) 

where G^g and Gg are the Green's function without and with 
the barrier present, respectively. Thus, in the absence of 
ambient flow and in the Fresnel number approximation, one has 
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* 


» 


IL » -10 log 10 |(l/2)(R/L) 2 [f 2 ([2N] 1/2 ) + G 2 ( [ 2N] 1/2 )j| 


( 20 ) 


Furthermore, for circumstances allowing the Fresnel number 
approximation, it is a good approximation also to set the 
factor R/L « 1, so one has 


IL • -10 log 10 j(l/2)[f 2 ([2N) 1/2 ) + g 2 ([2N] 1/2 )]j 


( 21 ) 


20 log 10 2[l + (2N) 1/2 ] 


N « 1 


s 20 log 10 2 + [20/tn 10](2N) 1/2 N << 1 

■ 6 + (12.3)N 1/2 dB N << 1 (22) 


-*• 10 log 10 (4ir 2 N) * 16 + 10 log 1Q N N >> 1 (23) 


which is a monotomically increasing function of Fresnel num- 
ber only. 

According to the analysis of the preceding section the 
magnitude of a Green's function when ambient flow is present 
is that of the Green's function when ambient flow is not 
present providing one lets x-*x/e, uj + uj/S, x g * Xg/s in the 
arguments of the latter. (This is true regardless of whether 
or not the barrier is present.) Consequently, the above 
approximate expressions for the insertion loss will still 
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apply to the case when there is an ambient flow, providing the 
Fresnel number is similarly transformed, i.e., 


N Cx L .yL,z L |xs,ys,ZsM 


N(x L /e,y L ,z L | x s/3,ys,Zs|u»/B) 


(24) 


In general, the ratio of the transformed and untransformed 
Fresnel numbers is spatially dependent. However, for the 
case when |y L /x L |«l and |y s /x s |«l, one has 


r L " 


|y L + \^y L \ 

l y s + x s /2y s i 
(r L + r s ) 2 = (y L + x 2 / 2 y L -y s - x 2 / 2 y s ) 2 


r s 8 


(25a) 

(25b) 



* ovn) 2 

• (y L • r s ) [x|/y s - x 2 /y L ] (26) 

L « 

[(t L - z s ) 2 » (y L -y s ) 2 

1 1/2 ( n - y s H*|/ys - *J/y L ] 


' 2 

[c 2 l' z s) 2+ (y L ' y s) 2 ] 1/2 




(27a) 

R - 1 

,(*L - *S> 2 * In - y s > 2 ] 


(x s *x L ) 2 

1 

2 | 

(z L -*s) 2 .(y L -y s > 2 j 1/2 


(27b) 


and, consequently. 


L _ R _ -(nAs^s ~ ( ys /y L )x L * 2x s*i ^ 

2 [(*L - Z S 3 2 * - 2 ] 1/2 

(Here, it should be recalled that the source and listener are 
presumed to be on opposite sides of the barrier, so y^ and yg 
have opposite signs. Consequently, the above gives L-R>0, 
as must be the case.) 

The above expression for L - R is bilinear in x g and x L 
so with the substitutions Xg + Xg/8, yg * one ^ as 
L - R -*• (L - R)/8 . Also the substitution w + w/B causes 
a/ 2 + 81/2. Consequently, in the case described above 

N -*■ N/8 3 

2 1/2 

Since 0 * (1 - M ) ' is less than 1, the transformed Fresnel 
number is larger than that corresponding to no flow. The 
insertion loss with ambient flow present is then given by 
Eqs. (21,22,23) only with N replaced by N/8 3 so one has in 
particular 


IL * 

6 ♦ (12.3)N 1/2 /B 3/2 

N<<1 

(28a) 

IL * 

10 log 10 (4* 2 N/8 3 ) 

N>>1 

(28b) 


In summary, the insertion loss is increased when there 


is an ambient flow, the increase being independent 
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of the direction of the flow. For larger values of the 
Fresnel number, the effect of ambient flow is to add an additional 
increment to the insertion loss of 


* 


% 




A(IL) = 10 log 10 [l/(l -M 2 ) 372 ] (29) 

If the Mach number is 0.5, for example, the additional 
insertion loss is 1.9 dB. 


G reen* s Function for Source Near Edge 

In Chapter 3, one of the limiting cases examined was 
that of a thin screen (v * 1/2) when ri‘g/L 2 -*-0, kr^r/L finite. 
This includes the case shown in Fig. 3 when the listener is 
many wavelengths from the diffracting edge and when the source 
is much closer to the edge than is the listener. However, the 
source is not presumed to be either very close or very far 
from the edge relative to a wavelength. The Green's function 
for this case can be constructed easily by the principle of 
reciprocity from Sommerfeld's known exact solution for the 
diffraction of an incident plane wave by a thin screen. (This 
was pointed out to one of the authors by Donald Lansing.) 

The result, for the case when the listener is in the shadow 
zone, is that the Green's function is given by 


G 


e ikL 

L 


nr I 


[f (X) - ig(X) ] 

t 
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(30) 

with 

X = [8rr s /AL] 1/2 |cos(?/2)| (31) 

The functions f(X) and g(X) are the auxiliary Fresnel functions 
described in the previous section. This Green's function may 
also be modified to take into account the presence of ambient 
flow by use of the transformation described previously. 

The fact that the above Green's function covers cases 
when rg is very close and not so close to the edge, that it 
is easily computed, and that it may be easily extended to 
include ambient flow suggests that it may be useful in studies 
of the diffraction of engine noise around wings while an air- 
plane is in flight. 

Sound from Distributed Sources 

In order to illustrate the application of the general 
theory to sound diffraction from a distributed source (Fig. 4), 
one may take, for simplicity, all the sources to be along the 
far side of the screen with 6g * 2* in each case and to each 
have Zg ■ 0 (i.e., the sources lie along a line transverse to 
the edge of the screen). The listener is considered also to 
have z coordinates equal to 0 and we consider r>>rg such that 


+ [f (X) - ig(X) ] 

c * 0 + e s 


SOURCE 



LISTENER 


Fig. 3. Limiting case of listener many wavelengths from 
edge of thin screen. 
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Fig. 4. Sketch of concepts utilized in the analysis 
of diffraction from a distributed source. 
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« 


v 


we may approximate L by r except in the exponent where it is 
taken as r+rg. In this manner, G reduces to 


G 


ST £ i l kr ** /4 I 

r 


jf(8r s /A] 1/2 cos [e/2]) 


- ig([8r s /A] 1/2 cos [e/2]) j e 


ikr, 


(32) 


If the source strength per unit length is taken as 
s(r s ) for a given frequency component, then the corresponding 
complex pressure amplitude in the far field is given by super- 
position as 

i [kr + ir/4] [ ikr Q 

P * /T f [f-ig]s(r s ) e b dr s (33) 

Here, for simplicity, we assume the source does not extend 
beyond the "trailing" edge, so all of the received sound is 
diffracted. 

For simplicity, we also assume the spatial extent of the 

source is somewhat less than a wavelength of the radiated 

ikr s 

sound such that we may approximate f, g and e by 
appropriate truncated power series expansions. If we do so, 
we have 


P 



1/2 - (4r g /A) 


1/2 


cos(0/2)e' in/4 ] 


t 
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• s(r s ) [1 + ikr g - (l/2)k 2 r|] dr s 


(34) 


One of the interesting aspects of the above is that the 
diffraction could enhance the received sound at low frequencies ♦ 
Suppose, for example, that the source were a quadrupole (e.g., 
as for jet noise) such that 



(35) 


Then the expression for p would reduce to 


P * -(l/4)k 2 | _/"rgS(r s )dr s 


- i lkr cos (6/2) e- 1 '/ 4 fr\ n 


r s S ( r s) dr s 

(36) 


The first term is weakened at low frequencies by the presence 
of the k factor while the second has a factor of k ' which 
may be larger when the frequency is low. The first term, 
incidentally, is just the sound field expected in the absence 
of the screen. [Time limitations, unfortunately, have pre- 
cluded a more thorough investigation of the question of 
whether diffraction could actually enhance the sound of a 
distributed source which tends to radiate as a quadrupole in 
the absence of a barrier.] 
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